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1.  INTRODUCTION  AND  SIMIARY. 

Most  numerical  approximation  processes,  such  as  interpolation,  quadrature, 
finite  difference  approximation,  finite  element  methods,  and  so  on, are  based 
on  enact  relationships  that  polynomials  (and  less  frequently,  trigonometric 
polynomials)  satisfy.  The  simplicity  and  variety  of  these  identities  have 
made  it  possible  to  construct  a vast  number  of  different  procedures.  These 
procedures  generally  do  very  well  in  a region  where  the  function  to  be 
approximated  is  analytic,  and  very  poorly  in  a neighborhood  of  a singularity 
of  the  function. 

This  paper  deals  with  the  approximation  of  analytic  functions  f on  an 
interval,  or  on  a contour.  These  functions  f may  or  may  not  have  a sin- 
gularity, i.e.,  a point  at  which  f'  does  not  exist,  at  the  end-points  of 
the  intervals  or  contours.  All  of  the  approximations  are  derived  via  the 
use  of  exact  relationships  that  the  function  C(f,h)  , i.e.,  Whittaker's 
cardinal  function, satisfies.  Corresponding  to  a function  f defined  on 
the  real  line  R , the  function  C(f,h)  is  defined  by 

00 

(1.1)  C(f,h)  - l f(kh)S(k,h) 

k*-<* 

whenever  this  series  converges,  where  h > 0 is  the  step-size,  and  where 

sin[?-  (x-kh)] 

(1.2)  S(k,h)  (x) -S 

£■  (x-kh) 


The  numerical  approximation  procedures  reported  in  this  paper,  and  which 
are  obtainable  via  the  use  of  S(k,h)0<J>  as  basis  functions,  where  4> 


denotes  a suitable  transformation  of  an  interval  onto  R , have  roughly 
the  same  accuracy  whether  or  not  the  function  to  be  approximated  has  a sin- 
gularity at  an  end-point  of  an  interval.  In  the  absence  of  singularities, 
this  accuracy  is  usually  not  as  good  as  that  obtainable  via  polynomial 
methods,  but  if  singularities  are  present,  this  accuracy  is  much  better 
than  that  of  polynomial  methods. 

The  function  C(f,h)  was  discovered  by  E.  T.  Whittaker  [55]  who  studied 
the  mathematical  properties  of  this  function  and  who  used  it  as  a means 
of  obtaining  alternate  expressions  of  entire  functions.  He  called  C(f,h) 

"a  function  of  royal  blood  in  the  family  of  entire  functions,  whose  distin- 
guished properties  separate  it  from  its  bourgeois  brethren".  The  study  of 
this  function  was  later  continued  and  considerably  extended  by  J.  M.  Whit- 
taker [56,57].  V-i  function  C(f,h)  then  played  an  important  role  in 
engineering  applications  in  the  transmission  of  information  as  a convenient 
approximation  of  f (Hartley  [16],  Nyquist  [30],  Shannon  [37]).  Engineers 
have  since  referred  to  C(f,h)  as  the  "band  limited"  or  "sine  function" 
expansion  of  f . 

The  mathematical  study  of  the  accurate  trapezoidal  formula  approximation 
of  the  integral  of  a function  f over  R , 


(1.3) 


00 

f(x)dx  ah  l f(kh)  , 
JR  k*-« 


developed  independently  of  the  study  of  C(f,h),  although  this  approximation  is 
identical  to  the  integral  of  C(f,h)  over  R.  It  was  Goodwin  [14]  who  seemed  to  be 
the  first  to  note  the  incredible  accuracy  of  this  formula  for  approximating 
the  integral  of  certain  functions  that  are  analytic  in  a strip  about  the  real 
axis.  The  incredible  accuracy  of  C(f,h)  and  of  the  trapezoidal  rule  as 
approximation  tools  in  the  family  of  functions  that  are  analytic  in  a strip 
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about  the  real  line  was  later  demonstrated  in  [25,26,27,54]. 

The  application  of  the  trapezoidal  formula  for  approximating  an  integral 
over  an  interval  other  than  R was  investigated  in  [28,34,39],  via  the 
use  of  transformations,  and  it  was  shown  in  [45,48]  that  the  most  effective 
transformations  are  those  that  are  a conformal  map  of  the  domain  of  analyticity 
of  the  integrand  onto  a strip  about  the  real  axis.  Later  [49],  this  trans- 
formation idea  was  used  to  construct  interpolation  and  approximation  formulas 
for  other  intervals  as  well  as  for  carrying  out  the  approximate  solution  of 
differential  [53]  and  integral  equations  [18,19,32,36,46,47].  Most  recently 

[52]  all  of  these  approximations  were  shown  to  have  the  optimal  functional 
-cn1/2 

form  0(e  ) for  the  rate  of  convergence  of  the  error  of  an  n-point 

approximation,  whether  or  not  the  function  being  approximated  has  singularities 
at  the  end-points  of  the  interval  (or  contour)  of  approximation. 

Although  the  present  paper  is  mainly  a summary  paper,  some  of  the  results 
in  it,  such  as  the  results  pertaining  to  the  case  where  a function  has  a 
singularity  on  the  interval  of  approximation,  are  new. 

The  function  C(f,h)  is  replete  with  beautiful  properties  and  formulas. 

The  known  properties  relevant  to  approximations  are  summarized  in  Sec.  2 of 
this  paper.  These  properties  are  basic  to  the  approximation  procedures  in 
later  sections  of  the  paper. 

In  Sec.  3 of  the  paper  we  define  a space  of  functions  that  are  analytic 
in  a strip  about  the  real  line.  In  this  space  the  identities  of  the  previous 
sections  are  no  longer  exact,  but  highly  accurate,  as  shown  by  the  error 
bounds. 

In  Sec.  4 some  of  the  approximations  of  Sec.  3 are  extended  to  an  arbitrary 
contour  by  use  of  conformal  mapping  [23,49].  Special  attention  is  given  to 
the  important  intervals  [0,1]  , [-1,1],  and  [0,®]  , and  examples  are 
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given,  of  approximation  rules  for  these  intervals. 

Section  5,  6,  7 and  8 deal  with  applications  of  the  previously  constructed 
approximation  procedures. 

In  Sec.  5 we  consider  the  approximations  of  important  transforms  over  the 
interval  (0,°°):  the  Laplace,  the  semi- infinite  Fourier,  the  Mellin  and 
the  Hankel  transforms  [22]. 

Sec.  6 deals  with  approximate  procedures  for  solving  ordinary  and  partial 
differential  equation  boundary  value  problems  [53].  It  is  here  that  the 
approximation  procedures  of  this  paper  are  particularly  powerful,  especially 
in  the  cases  where  the  singular  behavior  of  the  solution  on  the  boundary  is 
not  known  explicitly.  The  functions  S(k,h)*'t>  (or  product  of  these  for  the 
case  of  more  than  one  dimension)  are  the  basis  functions.  These  make  it 
possible  to  explicitly  write  down  highly  accurate  expressions  of  the  inner 
products  in  the  Galerkin  scheme  which  reduces  the  differential  equation 
problem  to  an  algebraic  problem,  even  for  the  case  of  nonlinear  differential 
equations.  Examples  of  the  approximate  solution  of  "model"  problems  illustrate 
the  application  of  the  method. 

In  Sec.  7 of  the  paper  we  apply  some  of  the  approximations  in  earlier 
parts  of  the  paper  to  the  approximate  solution  of  (singular)  integral 
equations.  Examples  are  given,  illustrating  the  approximation  procedures. 

Here,  too,  the  functions  S(k,h)0<J>  are  very  well  suited  for  easily  obtaining 
accurate  approximate  solutions. 

In  Sec.  8 we  summarize  the  main  ideas  used  for  the  implementation  of  the 
methods  on  a computer,  and  we  list  already  existing  computer  algorithms.  In 
addition  we  caution  the  user  against  possible  computational  pitfalls  resulting 
from  inaccurate  (or  improper)  numerical  evaluation  of  a function  in  the 
neighborhood  of  a singularity. 
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Sec.  9 discusses  rates  of  convergence  of  the  methods  of  this  paper,  and 

compares  these  with  rates  of  polynomial  methods.  It  shows,  moreover,  that 
-cn1/2 

the  0(e  ) rate  of  convergence  of  the  methods  of  this  paper  cannot  be 

improved  by  any  other  methods  of  approximation. 


2.  PROPERTIES  OF  THE  CARDINAL  FUNCTION. 


Just  as  polynomials  satisfy  certain  exact  relationships  in  the  space  of 
polynomials  of  the  same  degree,  the  same  is  true  of  the  cardinal  function  in 
a certain  space  of  entire  functions. 


DEFINITION  2.1.:  Let  h > 0 , and  let  B(h)  denote  the  family  of  functions 
f that  are  analytic  in  the  entire  complex  plane  (C  , such  that 

(2.1)  |f (z) | < Ce^ I z I /h  , 

2 

and  such  that  f e L (R)  , where  R = (-00,00)  . Let  k be  an  integer,  and 
let  us  set 

sin[£  (z-kh)] 

(2.2)  S(k,h)(z)  = . 

I (z-kh) 

If  f is  defined  on  R , the  Whittaker  cardinal  function  for  f with 
step-size  h is  defined  by 

oo 

(2.3)  C(f,h) (z)  = l f (kh)S(k,h) (z)  . 

k=-°° 


Set 

(2.4) 


jW 

Jk 


s(n)(j,D(k) 


n 


(£)  S(j  ,1)  (x) 


x=k 


In  particular 
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ff  • l 


1 if  j=k 
0 if  j*k 


(2.4a) 


0 if  j=k 


) / J 


if  jVk 


f /3  if  j=k 

4jy-  { .2(.1)k-j 

-A  ? if  j*k  . 

Oc-j) 

The  following  theorem  summarizes  the  known  explicit  relations  involving 
C(f,h)  and  S(k,h)  , in  the  case  when  f e B(h)  . 

THEOREM  2.1:  Let  f 6 B(h)  . Then: 


(2.5) 


f(z)  = C(f,h)(z)  for  all  z e I ; 


(2.6) 


J f(z)dz  = h l f(k,h)  ; 
Jr  v=-oo 


(2.7) 


|f(x)|2dx  = h l |f(kh)|2  , 

R W-oo 


and  the  set  (h  S(k,h)  )jc_.00  is  therefore  a complete  orthonormal  sequence 
in  B(h)  ; 

(d)  There  exists  a unique  function  g e L2  (-ir/h.Tr/h)  , such  that 


(2.8) 


2*  L 


e lztg(t)dt  ; 


(e)  The  function  g in  (d)  is  given  by 


, . - • . ■ * 


v 


(2.9) 


elxtf (t)dt 


r h l £(kh)elkhx 


. r tt  it 

if  - 7-  < X < - 

n n 


0 if  x > n/h  or  if  x < - ; 

h 


(2.10) 


i sin  £ (z-t) 

f(z)  = r f(t)  2 dt  ; 

TT  (Z't) 


(2.11) 


f'  e B(h)  ; 


(2.12) 


f(n)(kh)=h'n  l 5^  f(jh) 

j 3 -OO  J 


where  6^  is  defined  in  (2.4),  and  therefore,  by  (2.12), 

(2.13)  f(n)(x)=h*n  J ,(n)  f(jh)]  S(k,h)(x) 

k=-00  j = -00  j 


(i)  Let  g be  defined  as  in  (2.8).  Then 


(2.14) 


fCWt  - s 


and  in  particular f 


rx  fir  sin[(£  - k)C] 

(2.15)  S(k,h)(t)dt  - h[cr.  ♦ i 5-_ d£  ] 

o K * In  ^ 


where  (see  Table  2.1) 


(2.16) 


1 ("  sin  kg 

"In  ? ^ ’ 


fx 

Moreover,  if  f(t)dt  - 0 , and  if  f(t)dt  is  in  B(h),  then 

Jr  J -co 


" ,~.v  - 


-9- 


(2.17) 


00  00 


f(t)dt  = h l [ l Q .f(jh)]S(k,h)(x) 

k=-°o  i = -<x>  ■> 


(j)  Let  Pf  and  Pf  be  defined  by 


(2.18) 


(2.19) 


(Pf ) (x)  = lim 


(#f)(x)  = 


y-K> 


P.V. 

TTl 


m. 


+ 2-rri  JR  t-x-iy 


dt 


R 


f(t) 

t-x 


dt  . 


Then 


1 “ iir (x-kh)/h  , 

(2.20)  (PfHx)  - ^ J f(lth)  [ ■ 1 


and,  since  Pf  e B(h)  , 


(2.21)  G*)(x)  - I (j£(kh)  + 2il  l f (Jh)  [~ > S(k,h)(x) 
k=-°°  L * 1 j='oo  K J 


Similarly 


(2.22)  (Pf ) (x)  = 


00 

i l f(kh)  ~ (x-kh)S2(0,l)a[^(x-kh)] 
k-  * ® 


and  since  Hf  € B(h)  , 


00  00 


k-j 


(2.23)  (flfMx)  - i l ( l f (jh)  t-1'^ ]}S(k,h)(x) 

k=-oo  j = -00  J 


(k)  Let  a > 0 , and  let  a > - j . Then 


-10- 


f(t)dt 

■R  [(x-t)2+a2]a+1/2 

(2-24)  r?^~a  rTr/h  ■ . 

* T/P,  1.  J e'mg(t)|t|°K  (a|t|)  dt 

ir  T(a  + y)  ■’-tr/h  a 

where  g is  defined  in  terms  of  f , by  (2.8),  where  denotes  the 
Bessel  function. 


-1  r iair/2 . 


(2.25)  Ka(x)  = | Tr[sin(Tra)]"A[eJ-un/4J_ci(ix)-e'17ra/2Ja(ix)] 


and  where 


(2.26) 


00  r i -.m,  ,,.2ra+a 

J (x)  - l W (x/2) 

m=0  m!  F(m+a+l) 


In  particular 

(2.27)  Tk(a,a,h;x)  = 


S(k,h) (t)dt 

P 7 7 a+l/2 

R [(x-t)2+a2] 


172  f cos  [ £ - k) t ] t\  (aht)  It 

ir ' I'(a+l/2)  h h a 


and 


(2.28) 


Tk(a,a,h;&h)  =■  tk_£  (a,a,h) 


where 


(2.29) 


t,  (a,o,h)  * t2  I2ah) [ cos(kt)taK  (aht)dt 

tt ' rfa+1/2)  J0  a 


Hence 
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(2.30) 


* ~f~1vi/2  ' J.  fW>)Tlt(a,a.h;x) 

R [(x-t)2+a2]  k"'“ 


and  if  this  function  is  in  B(h)  , then 


(2.31) 


f(t)dt 


00  00 


, , 5*1/2  ‘J  <J  f(jh)T  (a>CI,h))S(k,h)(x) 
R t(x-t)2.a2l  k*-“ 


In  the  cases  when  a+1/2  >0  is  an  integer,  the  functions  x^(a,a,h;x)  can 
be  evaluated  explicitly.  In  particular,  if  a=l/2  , then  (2.30)  becomes 


t 

IT 


f(t)dt^ 
R [(x-t)  V] 


(2.32) 


- I f(kh)  ( y(1'e  Try/'h'lcos[TT(x-kh)/h]  » (x-kh)e  ^^sinfTrCx-khj/h])  . 
77  k=-°°  k (x-kh)  +y  ' 


(1) 


(2.33) 


2tt 


f(t)log[(x-t)2+y2]  dt 


oo  , 2 2 

£ l f (khK-y-iog  ir/h  ♦ [1-cos  (tt  (x-kh) /h}]  log  [I*:M  IT-  ] 

277  k—  l (x-kh)2 

♦ j*  l-(2-e  yt/^tl)cos{  (x-kh)t/h)  ^ j 


where  y denotes  Euler's  constant; 

(m)  Let  0 < a < 1 , and  let  g be  defined  as  in  (2.8).  Then 

r r/h 

|x-tr-*£(t)dt  - ■i.wvya^,,iw 

;R  2lT  J -ir/h 


(2.34)  f Ix-tl^fCtHt  • f/h  |tracct)e-ixtdt  , 

'R  ^ J-T r/h 
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and  in  particular 


(2.35) 


|x-t  |a  lS(k,h)  (t)dt  = harfa)cosCW2) 


t acos [(x-kh)t/h]  dt  , 


so  that 


(2.36) 


f Ix-tT'VtJdt  - h“ 

Jr  LK 


fir 


• l f(kh)  I t °cos [(x-kh)t/h]  dt 
v=  -co  J n 


These  results  follow  from  the  results  in  (k)  above; 


(n)  Let  g be  defined  as  in  (2.8).  Then 


(2.37) 


I log | x- t | f (t)dt 

Jr 


1 f 

4 


n/h 

-7r/h 


— 1 yf 

[gW-e  g(t)]  dt  _ 1 [Y+log(lT/h)]  g(0) 


and  in  particular 

f log  | x-t  | S(k,h)(t)dt 

Jr 

*-  jh  | Y+log(Tr/h) 


(2.38) 


r it 


l-cos[(x-kh)t/h]  dt  J 


so  that 


log|x-t|  f (t)dt 


(2.39) 


- ih  I f(kh)  ( Y+logO/h)  - r ?:cos[(x-M;)t/hl  dt  1 

z k— « <•  Jo  Z J 


This  result  is  obtained  from  (2.33),  by  letting  y •+  0 there. 
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TABLE  2.1 


simrx 

7TX 


Integrals  of  the  Sine  Function  from  0 to  n : a = | 

n J 0 


n 


i n t eg  r a l 


n 


dx 


integral 


1 

0.5394398722360836366+00 

51 

2 

0.451411 6667901403146+00 

52 

3 

0.533 09 323761 8271 9 84 6 +00 

53 

4 

0.474 96 96698836550 78 d+00 

54 

5 

0.520107164191 3085186+00 

55 

6 

U. 43320521 74977471 336+00 

56 

7 

0.51441 59971233052526+00 

57 

8 

0.43737422 50 5 78 199 73 6 +00 

58 

9 

0.5112301526369974586+00 

59 

1 0 

0.4898881711 538736606+00 

60 

1 1 

0.50919 574200821 66176+00 

61 

1 2 

0.491 56835 1663600880d *00 

62 

1 3 

0. 50773465781 25660746+00 

63 

1 4 

0.4927702093748031 35a +00 

64 

1 5 

0.5067436944720115796+00 

65 

16 

0. 4936724 1 51 782209536+00 

66 

1 7 

0.5059559079171768326+00 

67 

1 3 

0.494374552333662521 6+00 

63 

1 9 

0.5053297104401583626+00 

69 

20 

0.4949364995706954726+00 

70 

21 

0.5048226073041718776+00 

71 

22 

0.49539641 50843571536+00 

72 

23 

0.504 40 358519 80 504346 +00 

73 

24 

0.4957797661366431666+00 

74 

25 

0.5040515358439446466+00 

75 

26. 

0.4961041974385153476+0 0 

76 

27 

0.5037515950319091626+00 

77 

23 

0.4963823201653520646+00 

78 

29 

0.503 49 299327 7940413d +00 

79 

30 

0.4966233866311 80065d +00 

30 

31 

0.5032677369477718226+00 

31 

32 

0.4963343388552227596+00 

32 

33 

0.5030697632021874676+00 

33 

34 

0.4970204870276670436+00 

34 

35 

0.50289441 25557013266+00 

3 5 

36 

0. 49 7 18  5962 3 35 964 7 50d +00 

36 

37 

0. 50 2 73800 5 382 4382 94d +00 

37 

38 

0.497 33 40 2692690 66 31d +00 

88 

39 

0.502 59763321 5927650d +00 

89 

4a 

0.4974672909775606786+00 

90 

41 

0.5024709506908214816+00 

91 

42 

0.4975878678043931086+00 

92 

43 

0.502356043524901 3 35 d+00 

93 

44 

0.497697486705937335d+00 

94 

45 

0.50  2-251  35  66  77313093d +00 

95 

46 

0. 49779 757639 10 22 3 57d +00 

96 

47 

0.5021 555722141229696+00 

97 

48 

0.49 788 9 32 7564668 397d +00 

98 

49  - 

0.502 06 760 48 27451 950 d+00 

99 

50 

0. 49 79 73 74 0503081 8 27d +00 

100 

0. 501 98653516577951 46  + 00 
0.4980516616563337956+G0 
0.501911 5825934450376+00 
0.49812331 21215729286+00 
0.501 842079980753 1560+00 
0.4981908100179335646+00 
0.5017774537933251 556+C0 
0.4932531882343345326+00 
0.50171 72032611 157096+00 
0.498311 40 862927560 06 +C0 
0. 501 66091 25833090536+00 
0.4983658734334104546+00 
0.5016081906639415736+00 
0.4984169343055851546+00 
0. 501 55371 26935 1 1 5426+00 
0. 4934649019474760'’ 0 d+OC 
0. 501 51 21832244933926+00 
0.4935100478749203316+00 
0.5014683604663265636+00 
0. 49855261 43645059566+00 
0. 501  427001  5722  5036  56  + 00 
0. 49859281 633431361 76+00 
0. 501 337 90 36526 99 94 0d +00 
0.4986308454725350816+00 
0-.  501350 90 04573 34 9306 +G0 
0. 4986663732935892596+00 
0. 501 31 581 4563370501 6+00 
0. 498701 0537234768246+00 
0.5012325049928004736+00 
0.4987335252983492066+00 
0.5012503401830419506+00 
0.4937644130407322436+00 
0.5012207012502377436+00 
0. 4937933300630402586+00 
0. 501 191 9804979525946+00 
0. 498821 87897664961 16+00 
0.501164580131481 8556+00 
0.4988486530372720676+00 
0.501138411 1454785206+00 
0. 4 9 887 4 237 2 3097 1 99  3d +00 
0.5011133923582398766+00 
0.4988987091500927256+00 
0. 501089449570580461 d+00 
0. 4989221397841 47 706d+00 
0.5010665148309343176+00 
0.4989445942075476356+00 
0.5010445257913607426+00 
0.4989661321830880416+00 
0.5010234251415907156+00 
0.4989368086930455036+00 


3.  APPROXIMATIONS  OVER  THE  REAL  LINE. 


Whereas  the  relationships  of  the  previous  section  are  exact,  each  of 
the  formulas  (2.4),  (2.6),  and  (2.9)  provides  a method  of  approximation 
for  the  case  when  f does  not  belong  to  the  class  B(h)  . We  thus  intro- 
duce another  class  of  functions  defined  on  R = (-<*>, °»)  , for  which  the 
approximations  referred  to  above  are  extremely  accurate.  At  the  outset 
we  investigate  the  error  of  the  approximations  for  the  case  when  the  complete 
set  of  points  {kh}^=  oo  are  used  in  the  approximations.  We  then  also 

investigate  the  error  of  approximation  when  only  the  finite  set  of  points 
N 

{kh}^^  is  used,  and  h is  chosen  judiciously. 

DEFINITION  3.1.:  Let  d > 0 , and  let  denote  the  domain 
(3.1)  = (z  € I:  |£m  z|  < d)  . 


Let  p :>  1 , and  let  B^(Z^)  denote  the  family  of  all  functions  f that 
are  analytic  in  , such  that 

rd 


(3.2) 


-d 


|f(x+iy)  | dy  ■+  0 as  x -►  + ® ; 


and  such  that  N^(f  ,£>^)  < « , where 

N (f,0d)  - lim{(  f |f  (x+iy)  |pdx) 
p a y+d't  JR 


1/P 

1/P 


f , T)  VP  ) 

* ( |f(x-iy)|pdx)  > . 

JR  J 


*An  exception  occurs,  of  course,  for  the  case  of  Bq.  (2.10). 


(3.3) 


-15- 


If  p=l  , we  shall  simply  write  B(Pj)  and  N(f,Pj)  instead  of  B^Cty 
and  respectively. 

Let  us  set 

sin[r  (x-kh)] 

S(k,h)(x)  = ^ 

J-  Cx-kh) 


C(f.h)  = l f (kh)S(k,h) 

k=-00 

(3.4)  N 

CjjCf.h)  = l f(kh)S(k,h) 

k=-N 


E(f,h)  = f - C(f,h) 


EN(f,h)  * f - C^Cf.h) 


The  most  effective  application  of  the  formulas  of  this  section  occur  for 
the  case  when  f 6 B(d)  , and  when 

(3.5)  |f(x)|  < Ce'a|x| 

for  all  x € F , where  C and  a are  positive  constants. 


FIGURE  3.1  The  Region  V ^ of  Eg.  (3.1). 
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3.1  Error  of  Approximation  by  C(f,h)  . 

The  error  of  approximation  of  a function  f in  Bp(Pd)  by  C(f,h)  may 
be  expressed  explicitly  as  an  integral  [27]. 

THEOREM  3.1:  Let  f e B (ty  . Then  for  all  x e R , 


E(f,h)(x)  = 


2iri 


'R{ 


f (t-id~) 


(t-x-id)sin[^  (t-id)] 


(3.6) 


f(t+id  ) 


(t-x+id)sin[^-  (t+id)] 


dt 


This  result  forms  the  starting  point  for  obtaining  the  error  of  many  different 
types  of  approximations.  In  particular  the  following  theorem  was  proved  in 
[46]. 


THEOREM  3.2:  (a)  If  f € B(t?d)  , then 


1/2, 


N(f,fl.) 

(3.7)  2ird  H E(f,h)l(jo  , 2(id)^l£(f,h)l2  < 


(b)  If  f € , then 


1/2 


(3.8)  » E(f »h)l  2 , 2(7rd)  ' 8 £(£)«„  < 


N2cf>y 

sinh(^)- 


1 . 1 


(c)  If  f 6 Bp(t?d)  . l<P<»,|+±-l  , 

"Actually,  the  bound  on  lECf.h)^  is  new.  It  is  obtained  via  a direct 
application  of  the  inequality  > <t>l  2 < I f 1 2 » where  4>(x)  - R f(x-t)g(t)dt  , 

and  where  f e L2GR),  g e lV)  , to  (3.6). 


-17- 


then 


(3.9)  «E(f,h)l  < 


TT(2d) 


I7p 


[ 


r(^)r(i)  1/q  n (f,p.) 


4> 


i 


sinh(-^-) 


We  remark  that  once  we  have  bounds  on  IE(f,h)l-  and  IE(f,h)l  . we 

£ v 7 J 00  9 

jet  the 
inequality 


can  get  the  bound  on  lE(f,h)ls  , for  any  s between  2 and  « , via  the 


(3.10) 


■ E(f ,h)l s < lE(f ,h)» |/s  I E(f ,h)l^'2/S 


COROLLARY  3.5  (49):  ' Let  the  condition  (a),  (b)  or  (c)  of  Theorem  3.2  be 

satisfied  and  let  f satisfy  (3.5).  Then  by  choosing  h = [ird/  (aN)  ] 1//2 

,1/2 


(3.11) 


I EN(f  ,h)«s  < 


where  depends  only  on  f,  d and  a , and  for  appropriate*  s , as  in 
Theorem  3.2. 


* If  bounds  on  both  ■ E (f ,h) ■ 2 and  IE(f,h)lao  are  given  in  Theorem  3.2, 
then  (3.11)  holds  for  all  s e [2,«]  . 
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3.2  Error  of  Quadrature  by  the  Trapezoidal  Rule. 
Let  f e B(V^)  , and  let  us  set 


(3.12) 


and 


nN(f ,h) 


EN(f  ,h)  (x)dx 


N 

f(x)dx  - h l 
R k=-N 


fOch) 


(3.13)  n(f,h)  = lim  nN(f,h) 

The  error  n(f,h)  may  be  obtained  by  integrating  (3.6)  using  residues.  This 
yields  [22,26,27,45]  the  following  theorem. 


THEOREM  3.4:  Let  f € B(fl^)  , and  let  n(f,h)  be  defined  in  (3.14). 


(a)  Then 


(3.14)  n(f,h) 


sin(i  ft* id) ] 


f[t-id')«~td*lt)T/h 
sin[|  (t-id)] 


Moreover , 

1 e-7rd/h 

(3.15)  |n(f,h)|<  i S N(f ,P.) 

4 sinh(ird/h)  a 

If  in  addition  f satisfies  (3.5),  then  by  taking  h * [2ird/(oN)  ]1//2 

(3.16)  lnN(f,h)|  <Cie*C2irdctN^/2  , 
where  depends  only  on  f,  d and  a . 
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3.3  Fourier  Transforms. 


In  this  section  we  shall  give  a bound  on  the  error  of  the  approximation  used 
in  the  Fast  Fourier  Transform  method  [49]. 

Let  E(f,h)  and  EN(f,h)  be  defined  as  in  (3.4),  let  x e R , and  let 
us  set 


(3.17) 


«N(f,h) (x) 


f (t)elxtdt 
R 


' f (t)eixtdt 

R 


eixt  EN(f,h)(t)dt 


h l f(jh)eijhx  , |x| 
j-N 


. |x|  >£■  , 


< 


ir 

K 


and 


(3.18)  6(f ,h)  = lim  6M(f,h)  , 

N . 

Eq.  (3.17)  tells  us  that  the  sum  h £ f(jh)e1^  is  not  to  be  used  to 

j-N 

approximate  [ e*xtf(t)dt  if  |x|  > ir/h  . By  replacing  f(t)  by  f(t)e^xt 

JR 

in  (3.15),  we  get 


THEOREM  3.5:  Let  f 6 B(t?d)  , let  x€R  , |x|  < u/h  , and  let  5(f,h)  be 
defined  as  in  (3.19).  Then 


«cf,h)(x)  - JU  { iOziQi 


-(d-iu)(^x) 


JR  1 sin[^(u+id)] 


f(u-id')e 


-(d+iu)(J-  -x) 


sin  ^(u- id)] 


} 


du 


(3.19) 
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and  so 

(3.20) 

COROLLARY  3.6  [49] : Let  the  conditions  of  Theorem  3.5  be  satisfied,  let 
f satisfy  (3.5)  on  1R  , and  let  h * (ffd/aN)1^2  . Then 

1/2 

(3.21)  |6N(f ,h)  (x)  | < Cie'(7rdaN)  , |x|  < Tr/h 
where  is  a constant  depending  only  on  f,  d and  a . 


|6(f,h) (x)| 


1 -d% 
2 • . riTd^  6 

smh^) 


l*|) 


m 


roximation  of  Derivatives. 


In  applications  it  is  often  desirable  to  approximate  both  a function 
as  well  as  some  of  its  derivatives.  These  derivative  approximations  are 
readily  obtainable  by  differentiating  (^(f.h)  . Bounds  on  the  error  of 
approximation  of  f ^ on  R by  (^(f  ,h)  ^ are  readily  obtained  by 
bounding  the  integral  [ 2 5] 

E(f,h) (n)  (x) 


(3.22) 


2irl  R j=0  (n-j) ! (t-x-id)j+* 


f(t-id') 

sin[^(t-id)l 


j=0  (n-j)  ! (t-x+id)^+1  sin[^-(t+id)  ] 

which  follows  from  Eq.  (3.6).  The  details  of  bounding  E(f,h)^  and 
EN(f,h)^  are  carried  out  in  [251.  We  state  some  of  these  results. 


THEOREM  3.7:  Let  n > 0 be  an  integer,  (a)  Let  f e B(Pd)  , and  let 
ffd/h  > 1 . Then* 


(3.23)  1 E(f ,h) 


(n).  . n!  e1'2  N<f’V 

* ^ ~ ' 4 " I A ' 


2(Trd)  ' n hr  ,1/4  sinh(ird/h) 
ir  d 


(3.24) 


fnl  n!  e N(f,Pd)  (Tr/h)1 

>E(f ,h)  n • < - 


sinh(ird/h) 


See  the  footnote  on  page  16  re.  the  B E(f ,h) B 2 - bound  which  is 


new. 
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(b)  Let  f e B2(Pd)  , and  let  ird/h  > 1 . Then 


(3.25) 


»E(f,h)  (n)«2 


< 


n!  e N2(f,Dd) 
2tt 


(ir/h)n 

sinh(nd/h) 


and 

(3.26) 


!E(f,h) 


< 


n!  e1/2  N2(f,Pd) 

2 (ird) 1/2  (1  - -^-j) 
aL  dr 


1/4 


(Tr/h)n 

sinh(iTd/h) 


COROLLARY  3.8  Q25  ] : Let  f satisfy  the  condition  (a)  or  (b)  of  Theorem  3.7, 
and  let  f satisfy  (3.5).  Then  by  choosing  h = [trd/(aN)  ] 


EN(f,h) 


(n). 


(3.27) 


Jf^-c^f.h)  (n\ 


n+1 


< C1  N 


(mdaN) 


1/2 


for  any  s in  [2,“]  , where  Cj  depends  only  on  f,  d,  a and  n . In 
particular,  with  defined  as  in  (2.4), 

N rc+1  \I2 

(3.28)  |f(n)(kh)-h'n  l f(jh)6^|  < CtN"7”  e‘(7rdaN) 

j = -N  JK  = 1 


COROLLARY  3.9*:  Let  the  condition  of  Corollary  3.8  be  satisfied,  and  on  JR  , 
let 

(3.29)  |f(n)(x)|  < C2e'a|x| 

Then  there  exists  a constant  C^  depending  only  on  f,  d,  a and  n , such 
that 


*This  result  is  believed  to  be  new. 


(3.30) 


| f (n)  (x) 


n+3 


N N 1/2 

- h l ( l f ( jh)  6 j" ) S Ck ,h)  (x)  | < C^ff^e' (TrdctN) 


k=-N  j=-N 


= 3 


The  approximations  (3.28)  and  (3.30)  are  useful  in  the  solution  of 


differential  equations. 
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3.5  The  Indefinite  Integral. 


Let  us  present  an  approximation  of 


(3.31) 


I(x)  = f(t)dt 


in  terms  of  the  values  f (kh)  of  f [21 ]. 


THEOREM  3.10:  Let  f e B(Pj)  > and  let  g € B(P^)  where 


(3.32) 


s(x)  * IW  ' '<■> 


and  where 


(3.33) 


0 < 6 < 


On  ]R  , let 


(3.34) 


|f(x)|  < C e' 


and  let 


(3.35) 


a = min(a' ,2B) 


Then  for  h * [ird/(aN)  ]1/2  , 


f(t)dt  - f(t)dt 

e0  +e  0X  Jr 


(3.36)  N f N 

- h l I 

k=*-N  H«- 


IN  ( IN 

l { l Vi 

-N  N 3 


Jf(jh)- 


(eSj^e-Bih)' 


} S(k,I 


f(t)dt]  S(k,h)(x) 


< CXN  e‘(irdaN)‘ 
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where  is  defined  in  ( 2.1(J,  and  where  is  a constant  depending  only  on 
f,  d and  a . 

In  applications  it  usually  suffices  to  take  6 = 1/2  or  1 , changing  d 
instead. 


■ • . ' V.  ' ' . . 


T 
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3.6  The  Hilbert  and  Related  Transforms. 


Given  f € Bp(P^)  , p > 1 


each  of  the  integrals:  the  projection  integral 


(3.37) 


(Pf)(x)  = lira 
y+0+ 


1 

2iri 


f Ul 1 

JR  f-x-iy 


dt 


and  the  Hilbert  transform  integral 


(3.38) 


(Hf) (x)  = 


P.V. 

TTl 


' f(t) 
R *-x 


dt 


clearly  exist.  These  integrals  are  important  in  applications  dealing  with  Fourier 
transforms  such  as  in  the  solution  of  integral  equations  of  convolution  type.  We 
shall  give  approximations  of  these  integrals  in  terms  of  finite  sums  of  (2.20), 
(2.21),  and  (2.22)  and  (2.23)  and  we  shall  give  bounds  on  the  error  of  these 
approximations . 


THEOREM  3.11  [46,49]:  Let  f e B^ffl^)  , 1 < p < 00  , and  let  E(f,h)  be  defined 
as  in  (3.6).  Then 


(3.39) 


P[E(f,h)](x)  « - - .IS JlLL^d) 

JR  Hx-t-id)  sin  [ (t- id)  ] 


feiTTx/h+eiTr(t-t-id)/h 


(x-t+id)  sin  [ ^ (t+id)] 


2-  ft+idl  1 > 


dt 


and 


(3.40) 


H[E(f,h)](x) 


1 f [cosOrx/h)+e'iff(t'id)/h]f(t-id) 
^ ^ l (t-x-id)  sin  [ir(t-id) /h] 


iir(t+id)/h. 


+ rcos(Trx/h)+e^^nf(t+id)  | dt 

(t-x+id)  sin  [ir(t+id)/h] 
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Moreover,  if  p=l  , then 

(3.41)  »P[E(f,h) 13^ 

(3.42)  ltf[E(f,h)l<o 

while  if  p=2  , then 

(3.43)  IPIEff.h)]!. 

(3.44)  IHtECf.h)]!. 

Furthermore,  if  11 E (f  ,h) 0 2 is 

(3.45)  II  P[E(f  ,h)  ] B 2 


1+  \ e"ird/h 

N(f,Pd) 

2ird 

sinh (ud/h) 

1 + e -*d/h 
2nd 

N(f,Pd) 

sinh(ird/h) 

1 . | e -"d/h 

N(f,Pd) 

2 (mi) 1/2 

sinh(-rrd/h) 

1 + e _1Td/h 

N(f,Pd) 

2 (mi) 1/2 

sinh(7rd/h) 

bounded,  then 

* «H[E(f,h)]«2  = IE(f,h)l2  . 


In  addition,  if  p is  either  1 or  2, if  f satisfies  (3.5)  on  1R  , and  if 
1/2 

h = Od/(aN)]  , then  there  exists  a constant  depending  only  on  f,  d 


and  a , such  that 


(3.46) 


(Pf)(x) 


■yr  1 f (WO 
Z1  k=-N 


giirfx-kh)^ 

Tr(x-kh)/h 


< C1N1/2e‘(7rdotN) 


1/2 


(3.47) 


»(Hf)(x)  - ikINf(kh)^(x-kh)S2(0,l)=(^]ls  < CjN^V^^172 


for  all  s e [2,®]  . 
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We  remark  that  if  f satisfies  (3.5),  then  so  do  Pf  and  H f . Conse- 
quently, we  have  (see  (2.21)  and  (2.23)). 
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3.7  Singularities  on  the  Interval  of  Approximation. 

In  this  section  we  consider  the  extension  of  the  results  (k) , (£) , (m) 
and  (n)  of  Theorem  2.1  to  approximation  in  B(Pd)  • These  results  are  particu- 
larly useful  in  the  solution  of  singular  integral  equations  in  more  than  one 
dimension.  These  results  are  belived  to  be  new. 


THEOREM  3.13:  Let  f e B(t>d)  , and  let  Tk(a,B,h;x)  be  defined  as  in  (2.27). 
Then  for  1 < p < 00  , 1/p  + 1/q  = 1 , and  p(B+-y)  > j 


(3.50) 

< 


f(t)  dt 
[(x-t)2+a2] 


FI72 


l 


k=-00 


f(kh)Tk(a,S,h;x) 


„i/2ar 2(6H/2) 


dvq 


1 1/P 

r(P(B+2)  - 1/2) 

r(P(B+  1/2)) 


r(^) 

"nfT 


N(f,Pd) 

sinh(ird/h) 


Moreover,  if  f satisfies  (3.5)  and  if  tk(a,B,h)  is  defined  as  in 
(2.29)  then  there  exists  a constant  depending  only  on  f,  d,  a 


that  if  h = [ird/(aN)]1/2 
f(t)  dt 


[(x-t) 2+a2] 


8+1/2 


N ( N 'l 

l j l f(jh)x,  (a,B,h)[  S(k,h)  (x) 

:=-N  ' j*-N  > 


(3.51) 


(2.28)  and 
and  a such 


In  particular,  if  a*0  and  - j < B < 0 , then  for  f as  in  (3.50)  , 

I f f(t) |x-t| '2^-1dt  - h'28  r(-2B?c°s0rg)  £ f(kh)  rt28cos[(x-kh)t/h]dt 

I J]R  k=*  -<»  ' 0 

(3.52)  , 

rc-8)r(j*s)  N(f,od) 


sinh(TTd/h) 
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and  if  f satisfies  (3.5),  then  for  h » [ird/(aN)  ] ^ 

|[  fWI-tr26'1  dt  - h'2Bn:-28)cosW)  f I f f(jh) 
J]R  71  k=-N  ' j = -N 


t26cos[(k-j)t]dt) 


(3.53) 


•S(k,h) < Cx  N3/2e-(TrdoN) 


1/2 


where  depends  only  on  f,  6,  a and  d . 

Equations  (2.32)  and  (2.33)  yield  solutions  to  Laplace's  equation  in  the 
upper  half  plane.  Let  u satisfy 


(3.54) 


and  either 


(3.55) 


or 


uxx  + uyy=°  >xeIR»y>0 


(3.56) 


lim  u(x,y) 
y-K)+ 


li» 

y-0+  ** 


f(x) 


f(x) 


THEOREM  3.14:  Let  u be  the  solution  of  the  Dirichlet  problem  (3.54),  (3.55), 
where  f satisfies  the  condition  of  Theorem  3.2.  Then 

l _p*iry/h'|(.rtC  r-r.-uw'*  /wi 


ju(-,y)-  !l  7 f(kh)  ) cos [tt ( ♦ -kh) /h] * ( • ~kh) e ~ y sin[irl«-kh)/hl 

17  k-—  ' («-kh)2+y2  > 

< «E(£,h)l 


(3.57) 


where  IE(f ,h)lp  is  defined  and  bounded  as  in  Theorem  3.2.  Moreover,  if  f also 
satisfies  Corollary  (3.3),  then  by  choosing  h and  s as  in  Corollary  3.3, 


- 
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(3.58) 


'u(sy)~  Z {?  f(3h)y(1~e  3 } s(k,h)i 

71  k=-N  N (k-j)W  > 


< C,  N3/2e-C«)1/2 
“ 1 


where  depends  only  on  £,  d and  a 


THEOREM  3.15:  Let  f e B(Pd)  » and  let 


(3.59)  M(a,  f,  VA)  = [ | f(t+id‘) | + | f (t-id‘) | ] log[ (t-ar+d‘]dt  < « 

a Jr 


for  all  finite  a?R  , Then  the  function  u=u(x,y)  , 


(3.60) 


u(x,y)  = j-  f log[(x-t)2+y2]f(t)dt 


which  solves  the  Neumann  problem  (3.54),  (3.56)  satisfies 


,°°  , mi  2 2 

Ju(*,y)~yr  l f (khM  -Y-logiT/h+[l-cos{Tr(*  -kh)/h}]  log[  - -Z  ] 

/7r  k— o»  C--kh)^ 


(3.61) 


♦ f”  l-(2-e'ty/h)co»[(--kh)t/linjt. 

. fl  t J 


M(a,  f,  0d)e"ird/h+4(h/d)N(f,0d) 
8Trsinh(ird/h) 


Moreover,  if  f satisfies  (3.5),  and  if 


(3.62) 


f(t)dt  - 0 , 


then  for  h ■ [Trd/(aN)  ' ] , 


j 

4* 


. -32- 

k=-N  lj-N  (k-j)2 

j^k 


(3.63) 


f l-(2-e~ty/hJcosUk-j)tldt  1 Sfk  h), 

Jo  t > 


< Cx  N3/2e‘(irdaN) 


1/2 


where  depends  only  on  f,  d and  a . In  particular,  if  y=0  , (3.63)  re- 
duces to 


u N , N 
log|.-t|f(t)dt-|-  l l f(jh) 
R k=-N  'j=-N 

j#k 


at  }s(k,h) ll 

0 t > 


(3.64) 


1/2 


4. 


FORMULAS  OVER  FINITE,  SEMI - INFINITE  INTERVALS  AND  CONTOURS. 


In  this  section  we  shall  consider  the  extension  of  some  of  the  formulas 
of  the  previous  section  to  finite  and  semi  infinite  intervals,  and  more  generally, 
to  contours  [45,49]. 

The  success  of  the  methods  of  the  previous  section  for  functions  in  V ^ 
suggests  that  a problem  over  an  arbitrary  interval  be  transformed  into  one  over 
(-oo,ao)  in  such  a way  that  the  conditions  of  the  theorems  of  the  previous  section 
are  satisfied.  In  this  way,  some  formulas  over  contours  other  than  (-“,»)  may 
be  obtained  directly.  However,  not  all  formulas  transform  directly,  and  we  must 
make  certain  simple  "adjustments"  in  order  to  make  direct  transformation  possible. 
The  following  definition  is  fundamental  for  the  remainder  of  the  paper. 

DEFINITION  4.1:  Let  V be  a simple  connected  domain,  with  boundary  dV  , let 
a and  b*a  be  points  of  dV  and  let  V ^ be  defined  as  in  (3.1).  Let  $ be 
a conformal  map  of  V onto  V d , such  that  <|>(a)  = -«  , <J>(b)  * 00  • Let 
ip  = <p  1 denote  the  inverse  map , and  set 

(4.1)  T * {ip (x)  : -oo  < x < ® } 

Given  4>  and  \p  , we  denote  by  z^  = z^(h)  the  points 

(4.2)  z^  a <J>(kh)  , k*0,±l,±2,. . . 

Let  B(P)  denote  the  family  of  all  functions  F that  are  analytic  in  V , such 
that 


(4.3) 


| F(z)dx|  -*-0  as  u -►  ±oo 

^(u+L) 
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wherc 

(4.4)  L = { iy  : y is  real  , |y|  < d}  , 
and  such  that 

(4.5)  N(F,P)  = lim  inf  |F(z)dz|  < °° 

cr30  cic°  Jci 

We  remark  that  if  Fe  B(0)  , then  f defined  by 

(4.6)  f = [F°^]^' 

is  in  8(0^)  as  defined  in  Definition  3.1. 

Let  us  next  give  four  commonly  used  transformations  <p  and  the  corresponding 
inverse  functions  ^ , corresponding  intervals  r = [0,1], [-1,1]  and  . 

EXAMPLE  4.1:  r = [0,1]  - In  this  case 

(4.7)  V = (z  : |arg  — | < d)  , 

1-z 


FIGURE  4.1  The  Region  V of  Ex.  4.1. 

The  boundary  of  V consists  of  two  circular  arcs  which  intersect  with  angle  2d 
at  0 and  at  1 . 
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The  function  p , the  inverse  function  and  the  points  z^,  are  given  by 
w = 4>(z)  = logj— <=»  z = y + ^ tanh^ 


(4.8) 


1 . 1 


zk  = ^ + J tanh(kh/2)  , k=0,±l,±2,. 


EXAMPLE  4.2:  T = [-1,1]  . In  this  case 

(4.9)  V = (z  : |arg(~|)|  < d)  . 


The  functions  <J>  and  ip  and  the  points  Zj.  are  given  by 
w = 4>(z)  = log(p|)<=»  z = i|»(w)  = tanh ~ 


(4.10) 


z^  = tanh(kh/2)  , k=0,±l,±2,. . . 


FIGURE  4.2  The  Region  V of  Ex.  4.2. 

We  shall  give  two  examples  for  the  case  T = [0,°°]  . The  first  of  these  is 
useful  in  the  case  when  the  function  f to  be  approximated  is  analytic  in  a 
sector,  |arg  z|  < d , while  the  second  is  useful  if  f is  analytic  only  in  a 
strip  of  width  2d  symmetric  about  the  real  axis  (more  precisely,  in  the  region 
V of  Ex.  4- -see  Figure  4.4). 


EXAMPLE  4.3: 


r = [OH  . In  this  case  V is  the  sector 


(4.11) 


V = (z  : |arg  z|  < d } 


FIGURE  4.3  The  Region  V of  Ex.  4.3. 


The  functions  4>,\p  and  the  points  z^  are  given  by 


(4.12) 


w = ip(z)  = log  z «=»  z = ip(w)  = e 
lcVi 

Zj^  = e"1  , k=0,±l,±2,. . . 


w 


EXAMPLE  4.4  [22] : r = [OH  . In  this  case 


(4.13) 


V = (z  : |arg  sinh(z)  | < d]  , 0 < d < ir/2 


FIGURE  4.4  The  Region  V of  Ex.  4.4. 
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The  functions  <p,\p  and  the  points  z^  are  given  by 


w = $(z)  = logfsinh  z]  z = iHw)  = log[eW+/  l+e2w  ] 


,kh 


(4.14)  z,  = log[e  +/  1+e  ] , k=0,±l,±2,.. 


0kh  1 3kh  x 3 5kh  5 7kh  kh  . . 

-e  -5e  • me  < A 


kh 

The  expansion  in  (4.14)  is  preferred  if  e < .1  since  the  formula 


z^  = log[e^+/  l+e“  ],  while  mathematically  exact,  is  computationally  inaccurate 
for  small  e^1  . In  applications,  the  accurate  computation  of  f(z^)  » where 
f(zk)  is  to  be  approximated, is  important,  especially  near  a singularity  of 
f (see  Sec.  8.2) . 

We  shall  next  review  -the  known  results  corresponding  to  those  in  Sec.  3 
which  may  be  extended  to  approximation  over  a contour  r as  described  in  Def.  4.1, 
albeit  some  minor  modifications.  These  modifications  are  described  in  the 
theorems,  as  we  shall  present  them. 

The  condition  (3.5)  takes  on  a simple  general  form:  let  g be  defined  on 
T , and  let 

(4.15)  |g(x)|  < Ce"al*Wl 


for  all  x € r , where  C and  a are  positive  constants. 

Two  identities  play  an  important  role  in  obtaining  all  of  our  bounds.  These 
are  described  in  the  following  theorem. 


THEOREM  4.2  [49] : Let  F e B(D)  . Then  the  identity 
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(4.16) 


F(x) 
<p'  (x) 


7 F(zk) 

l S(k,h)o<j>(x) 

k=-®  <p'  (zk) 


sin[TT<ft(x)/h]  f F(z)dz 

2rri  ^ dV  [0  (z) -4>  (x)  ] sin  [tt4>  (z)  /h] 


is  valid  for  all  x € r . Moreover 


(4.17) 


00 

F(x)dx-h  £ 

F k=-°° 


eXp[iH^I  Sgn  ImcKz)] 
sin[n4)(z)/h] 


F(z)dz 


We  remark  that  (4.17)  is  obtained  from  (4.16)  by  multiplying  (4. 16)  by 


<p ' (x)  and  integrating  over  r . 
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4.1  Interpolation  over  r . 


THEOREM  4.3  [45,49] : Let  <fr’F  € B(t>)  . Then 

00 

(4.18)  |F(x)-  X F(zJS(k,h)«*(x)|  < — 

k— • 2ird  sinh(Trd/h) 

for  all  x € r . Moreover,  if  F is  bounded  on  V by  the  right-hand  side  of 
(4.15),  then  by  taking  h = [>d/(aN)  ]1//2 

(4.19)  |F(x)-  [ F(z.)S(k,h)o^(x)|  < 

k— N K “ 1 

for  all  x e r , where  C-^  depends  only  on  F,  d and  a . 

It  may  happen  that  F 6 B(D),  but  <j>'F  £ B(£>)  . If  the  limits 

lim  F(x)  = F(a) 
x-*-a 

(4.20) 

lim  F(x)  = F(b) 
x-*t> 


exists  and  are  bounded,  where  the  limits  are  taken  along  T , then  it  may  also 
be  the  case  that  <t>'G  e B(P)  , where 


(4.21) 


F - 


-*54> 


e*i<P 

M + e'1^ 


F(b) 


This  device  is  often  useful  in  applications. 

Let  us  next  illustrate  the  formula  (4.19)  and  (4.21)  for  the  case  of  the 
transformations  in  Examples  4.1  to  4.4. 
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EXAMPLE  4.5:  For  the  case  when  r = [0,1]  , let  F be  analytic  and  bounded 
in  the  domain  V of  (4.7).  On  [0,1]  , the  condition  (4.15)  becomes 


(4.22) 


|F(x)|  < C xa(l-x)a 


If  F satisfies  (4.22)  on  [0,1]  , then  by  taking  h = [ud/ (aN) ] ^ , (4.19) 
becomes 


F(x)-£sin{ 


1 x>  k=-N 


N (*l)KF(zk) 


(4.23) 


< c/ZV^ 


log 
1/2 


where  z^  is  given  by  (4.8),  and  where  depends  on  F,  d and  a . If  F 
does  not  vanish  on  0 or  at  1 , the  function  G of  (-4.21)  may,  provided  that 
the  limits  (4.20)  exist.  The  function  G takes  the  form 

(4.24)  G(x)  = F(x)-(l-x)F(0)-xF(l) 

We  remark  that  if  F e B(t0  and  if  F 6 Lip^flJ) , where  V denotes  the  closure  of 
V , then  (J>'G  <=  B(0)  and  G satisfies  (4.22)  on  [0,1]  [ 36  ] . The  formula 
(4.23)  does  a good  job  of  interpolating  functions  such  as 

F(x)  = x1Z^(i-x) *^Z2^0g  x * or  p(x)  - sin(Trx)log(l-x) , etc.  on  [0,1]  . 

The  case  of  T = [-1,1]  is  similar  to  that  of  T = [0,1]  . 


EXAMPLE  4.6:  Consider  the  case  of 
F be  analytic  in  the  sector  (4.11). 

(4.25)  ] F(x) | < { Cx“ 

( Cx  a 


V as 

in  Ex.  4.3, 

with  T * [0,®]  . Let 

On 

[0,«]  , let 

F satisfy 

if 

0 < x < i 

if 

1 < x < » 

1 

I 


a condition  which  is  equivalent  to  (4.15).  Then  by  taking  h = [7rd/(aN) , 


(4.19)  becomes 


(4.26) 


|F(x)-£  sin{J-  log  x}  l F(e — 2- 

k=-N  log  x-kh 


< C1N1^2e"^1T<^2" 


where  depends  only  on  F,  d and  a . If  F does  not  vanish  at  0 or  at 
00  and  if  the  limits  (4.20)  exist,  then  we  may  be  able  to  effectively  apply  (4.26) 
to  G (see  Eq.  (4.21))  where 


(4.27) 


G(x)  ■ F(x)^F(0)-^F(»)  . 


The  formula  (4.26)  does  an  accurate  job  of  interpolating  functions  F such  as 
F(x)  * x2/,3(log  x)/(l+x)  , or  for  F(x)  = x5/2e‘xsin  x/2  , etc. 


EXAMPLE  4.7:  Let  us  again  take  f = [0,°°]  , for  the  case  of  Ex.  4.4.  For  this 


case  the  condition  (4.15)  becomes 


(4.28) 


|F(x)|  < {C 


if  0 < x < 1 

ax  ic 

if  1 < X < °9 


If  ’ F e B(P)  and  if  F satisfies  (4.28)  on  [0,°°]  , then  by  taking 
h = [ird/ (oN)  ] l^2  , we  get 


h N C-1)KFCzJ 

F(x)-£  sin  {£-  log[sinh  x]}  7 

k»-N  log[sinh  x]-kh 


(4.29) 


< CXN1/2  e'CirdaN)‘ 


for  all  xe  [0,<»]  , where  is  defined  in  (4.14)  and  where  depends  only 


•n  ■/  -- 


♦ 
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on  F,  d and  a . If  F is  analytic  in  V but  if  F does  not  vanish  at  0 
or  at  ® , then  the  function  G may  satisfy  the  conditions  of  Theorem  4.3 
well  as  (4.28),  where 

(4.30)  G(x)  - F(x) F(0) ^L£f(«)  . 

1+sinh  x 1+sinh  x 


The  formula  (4.29)  does  well  at  interpolating  functions  over  [0,®]  which  may 
be  oscillatory  on  (0,®)  , but  which  may  have  a singularity  at  x=0  . For 
example,  (4.29)  does  well  at  interpolating  F(x)  = xalog[l- (-s^  x)2]  e'0^  , < 
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4.2  Quadrature  over  r . 


Eq.  (4.17)  yields  the  following  theorem. 


THEOREM  4.4  [45]:  Let  F € B(fl)  . Then 


(4.31) 


“ FCzk) 

F(x)dx-h  l — 

r k=-o°  <j>'  (zk) 


-7rd/h 


2 sinhfrrd/h) 


N(F,t?)  . 


Moreover,  if  F/<J> 1 satisfies  (4.15)  on  T , then  by  taking  h = [27rd/(aN)] 

,1/2 


1/2 


(4.32) 


N F(zJ 

F(x)dx-h  l 

k=-N  $’(zk) 


< C1  e 


- (2iTdaN) ' 


where  depends  only  on  F,  d and  a 


EXAMPLE  4.8:  Let  r = [-1,1]  , and  let  F e B(P)  where  V is  defined  in  (4.9) 
On  [-1,1]  , let 

->  ot-1 


(4.33) 


| F(x) | < C (1-x4)  , a > 0 , C.  > 0 , 


a condition  equivalent  to  F/<j>'  satisfying  (4.15)  on  F 
taking  h = [2nd/(oN)  ]1//2  , 

rl  N 


[-1,1]  . Then  by 


. kh  kh_, 

F(x)dx-h  l — 0 

-1  k=-N  (1+e  j L e^+l 


(4.34) 


< cx  e‘(2irdctN) 


1/2 


where  depends  only  on  F,  d and  a . The  formula  (4.15)  yields  accurate 
results  for  the  integration  of  function  F such  as  F(x)-(l-x) *^(l+x)  ^5log(l-x) 
or  F(x)-(l-x)*exp{-2/(l-x)}  . 
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EXAMPLE  4.9:  Let  T = [0,°°]  , and  let  F € B(fl)  where  V is  defined  in  (4.11) 
On  [0,-]  , let 


C x01’1 

(4.35)  |F(x)|  <{c  x-«-l 


, 0 < x < 1 

> 1 < x < °° 


a condition  which  is  equivalent  to  F/<J>'  satisfying  (4.15).  Then  by  taking 
h = [2iTd/(aN)]1/2 


(4.36) 


F(x)dx-h  l e^Fte^1) 
k=-N 


< ci  e 


-(2TTdoN)1/2 


where  depends  only  on  F,  d and  a . The  formula  (4.36)  does  an  accurate 
job  of  integrating  function  F such  as  F(x)  = xa"1/(l+x)2ct  , or 
F(x)  = x ^2sin(x/2)e  x . 

EXAMPLE  4.10  [22] : Let  r = [0,°°]  , and  let  F € B(P)  , where  V is  defined 
in  (4.14).  On  [0,°°]  , let 


(4.37) 


|F(x)|  < 


C xa'1  if  0 


C e'ax  if  1 < x < 


< x < 1 1 

" “ ) “ 


> 0 


a condition  which  is  equivalent  to  F/<p'  satisfying  (4.15).  Then  by  taking 
h - [2Trd/(oN)]1/2  , 


. 

F(x)dx-h  J F jlogrekh+  \/l+e2kh  ] } 

k-N  /~~=m  l V ) 


1+e 


(4.38) 


< Cx  e 


-(2irdoN)1/2 


where  depends  only  on  F , d and  a . The  formula  (4.38)  does  an  accurate 

job  of  integrating  functions  F such  as  F(x)  * x~1//2log[l-?—  -x)eX//2  , 

2 x 
F(x)  - x'2/7e*x  , or  F(x)  - x'5exp{-[(x-S)2+2]1/2-l/x2}sin(3x)  . 


I 
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4.3  Approximation  of  Derivatives  on  r . 

Except  on  ]R  , the  formula  (4.18)  is  not  useful  for  accurately  approximating 
derivatives  of  F on  r , since  the  terms  <t>'S(k,h)°<l>  are  unbounded  on  r . 

Then  to  get  a formula  for  approximating  f^  on  r it  becomes  necessary  to 
modify  (4.18)  by  introducing  a "nullifier"  function  g with  the  property  that 

|g(x)S(k,h)°  <j>(x)| 

is  bounded  on  r , for  n=0,l,...,m  . Upon  replacing  F in  (4.16)  by  F<j>'/g 

we  get 


THEOREM  4.5  [ 25] : 


(4.39) 


Let  F<J)'/g  € B(P)  . Then  for  all  x e r , 
" F^zk^ 

F(x)  - l £-g(x)S(k,h)o*(x) 

k— " gCzjj) 


= g(x)sin[ir<t>(x)/h] 
2tti 


F(z)<fr'(z)/g(z)  dz 
|>(z)-<J>(x)  ]sin[Tr<j>(z)/h] 


Let 


(4.40) 


f dj 

n 

g(x)  sinrir<Kx)/hl 

[Xx\ 

4>(z)-<j»(x) 

for  n=0,l,...,m  , for  all  x € r and  z 6 3P  , where  C 2 is  a constant 

depending  only  on  m,  g and  F . Then  there  exists  a constant  , depending 
only  on  m,  g,  d and  F , such  that  for  all  x e r , 


(4.41) 


F(n)  (x)- 


d 

■ar, 


n 


JL 


g(x)S(k,h)  °4>  (x) 


< C3h’n  e'^ 


f n-0 | 1 f • • • ,m 


Moreover,  if  F/g  satisfies  (4.15)  on  r , then  by  taking  h = [7rd/(aN)] 
then  for  all  x e r 


1/2 

9 


(4.42) 


F(n;)  (x) 


N F(z.) 

l g(x)S(k,h)o<j)(x) 

k*-N  g(z^) 


n+1 

< C, 


-(TrdctN) 


1/2 


n=0,l , . 


,m 


where  depends  only  on  m,  g d,  F and  a . 

The  function  g takes  on  different  forms  for  different  values  of  <J> 
although  g(x)  = [l/<)>'(x)]m  is  usually  satisfactory.  In  addition  Theorem  4.5 
presupposes  that  all  derivatives  of  F vanish  at  a certain  rate  as  x-*-a  and  as 
x+b  along  r . We  shall  present  special  forms  for  g as  well  as  procedures 
for  making  a number  of  derivatives  of  F vanish  at  a and  at  b in  the 
examples  which  follow. 


EXAMPLE  4.11  [25] • Let  r = [-1,1]  , and  let  V be  defined  as  in  (4.9).  Let 
us  take 


(4.43)  g(x)  = (l-x2)m 

and  let  us  assume  that  [F<J>'/g]  e B(P)  where  F(x)4>' (x)/g(x)  = 2F(x)(l-x2)  m’1' 
Furthermore,  on  [-1,1]  , let  us  assume  that 


(4.44) 


F(x) 

F(x) 

g(x) 

(l-x2)1" 

< C(l-x2)a 


1/2 


a condition  on  F/g  equivalent  to  (4.15).  Then  by  taking  h = [7rd/(aN)] 

fnl  (d)n(h  2 m 7T  1+x  y (-DkF(zk)/(l-zk)m  ) 

pCn)(x)-  4 ?d-x  ) sin[£log(^)]  l j— ^ *_} 

iaxJ  l * nix  log[£|]-kh  > 


< Cj  rrv^172 


p n-0 y 1 f • • • ,m 


(4.45) 
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for  all  x e r . More  generally,  let  be  analytic  and  bounded  on  V and 
F^  e Lip^CI?)  , where  V denotes  the  closure  of  P . In  this  case  (4.44)  may 
not  be  satisfied.  However,  the  function 


(4.46) 


G = 


F 


- P 


m 


satisfies  (4.44)  as  well  as  all  of  the  other  conditions  required  in  (4.45),  where 

p is  constructed  as  follows.  Set 
m 


(4.47) 


po(x)  = ao 


1-x 

2 


+ bo 


1+x 

2 


where 


(4.48) 
and 

(4.49)  Pk+1M 


aQ  = f(-l)  , bQ  = f(l) 


Pk(x)  + Wt 


fl-y)k+2 

- ' ' r-  -| 


fl+x 

ik+1 

l-x 

|k+!| 

1+x 

2 j 

+ Vl 

hn 

[~r 

k+2 


where 


(4.50) 


Vi  ■ 

2k+1 

[f(k+l)(-l) 

- pk(k+1)(-l)] 

(k+1) ! 

hfc*!  ‘ 

(-2)k+1 
Ck+1) ! 

[f(k+1)(-l) 

- Pk(k+1)(l)] 

EXAMPLE  4.12  [23] : Let  T * [0,«]  , and  let  F be  analytic  in  the  region  V 
of  Eq.  (4.11).  Let  us  assume  that  exists  on  [0,°°]  , let  us  take 

g(x)  * x"1  , and  let  us  assume  that  4> ' F/g  e B(D)  , where  [<t>'F/g](x)  - F(x)/xm+1 

Furthermore,  on  [0,®]  , let  us  assume  that 
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(4.51) 


F(x) 


C xa  , 0 < x < 1 


x° 


C x , 1 < x < oo  , n = 0,  1, . . . , m . 

a condition  on  F/g  equivalent  to  (4.15).  Then  by  choosing  h = [iTd/(aN)] 
I’d  'n 


1/2 


FCn) (x)- 


(4.52) 


3ST 


. t _ „ N ( T,kcr  kh.  -mkh 

t {Aintf  log  x]  I (-1)  ^ 

k=-N  log  x-kh 


ir 

n+1 


< Cx  N*1"  e't7rdaN)1/2  , n=0,l,...,m 


for  all  x s [0,«]  , where  Cj  depends  only  on  F,  d,  m and  a . More 
generally  if  F^  <=  B(D)  and  if  F^  e Lip^fP)  , then  F may  not  satisfy 
(4.51).  However,  the  function  G usually  does  satisfy  all  of  the  requirements 
of  (4.52),  where 


G(x)  = F(x)  - e X l a,  x1' 
k=0  K 

a0  = F(0) 


m 


(4.53) 


ak  - [F°°(0)  - l ^)k  ^ a ] , k-1,2 m 

K j=0  (k-j )!  J 


J ■"  • • * 
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4.4  Approximation  of  the  Indefinite  Integral  on  r . 


We  shall  give  a general  formula  for  approximating 


(4.54) 


100 


rx 

F(t)dt 

a 


x e r , 


in  Theorem  4.6  below.  The  special  forms  of  this  formula  for  the  cases  of 
Examples  4.1  to  4.4  are  omitted,  since  these  special  cases  are  simply  obtained 
by  direct  substitution.  The  results  of  this  formula  are  especially  suited  to 
the  solution  of  linear  initial  value  problems  and  to  some  linear  Volterra  integral 
equations..  The  solutions  of  these  equations  may  have  a singularity  or  a boundary 
layer  at  one  or  both  end-points  of  T . 


THEOREM  4.6  [24] : Let  F e B(P)  , and  along  T 

(4.55)  j —■(?■)■  j < c e~a'  1 

U’(x)  1 


let 


where  C and  a'  are  positive  constants.  Let  0 < 8 < v/d  , and  let  G € B(P)  , 
where 


(4.56) 


G(x)  = 


rX 

F(t)dt  - 
a 


e*tf$(x) 

e*5S<t>(x)+  ^341  (x) 


b 

F(t)dt 

a 


Let 

(4.57)  a = min(a' ,6) 


1/2 

f 


Then  for  h = [Trd/(aN)] 
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(4.58) 


WCx) 


F(t)dt  'j64>(x)  li^Cx) 

A T A 


F(t)dt 


N f N r r(.z.j  , - v 

J-N  iji-N  ^ L^rt  ' Te^Te^W  1 J } ' S(k'h,°*W 


F(z.) 


<C:  N e'(7TdaN} 


1/2 


for  all  x e r , where  ak  is  defined  in  (2.15),  and  where  depends  only 
on  F,  d and  a . 


In  applications  the  integral 


F(t)dt  is  approximated  by  means  of 


formula  (4.32),  i.e.,  by  h £k=_N  F(zk)/$’ (zk)  . 

For  example,  in  applications,  the  formula  (4.58)  may  be  used  to  approximate 


an  integral  such  as 


[x  ,-2/3 
0 


(log  t)  (l-t)'5//4dt  , x€  [0,1] 
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4.5  Singular  Integrals  on  r . 


In  this  section  we  consider  the  approximation  of  the  integrals 


(4.59) 


(ffF)(x)  = P.V. 


f £01 


t-x 


dt  , x € f ; 


(4.60)  (l/ftF)(x)  = 


[(x-t)2+y2]'6'1/2  F(t)dt  ,xer,y>0,B>  -1/2  , 


(4.61)  (PgF) (x)  = 


x-t|6'1  F(t)dt  , x e r , B > 0 ; 


(4.62) 


(WF)(x)  i 


log|x-t|  F(t)dt  , x € r . 


The  results  of  this  section  are  believed  to  be  new. 

Although  the  last  three  of  these  integrals  can  be  evaluated  via  the  method 
in  Sec.  4.2,  by  splitting  up  each  integral  as  an  integral  from  a to  x plus 
an  integral  from  x to  b , the  methods  developed  in  this  section  are  more 
efficient,  since  we  derive  explicit  approximate  expressions  that  are  valid  for 
all  x e r . This  increased  efficiency  is  especially  important  in  the  solution 
of  integral  equations,  where  the  major  difficulty  is  the  evaluation  of  many 
singular  inner  product  integrals. 

Although  the  case  of  ffF  is  done  fairly  generally,  due  to  its  importance, 
the  procedure  for  the  case  of  the  remaining  integrals  is  illustrated  only  for 
special  intervals.  We  have  also  left  out  the  case  of  pF(t)log[ |x-t|  +|y|“]  dt  , 
since  formulas  for  this  case  can  readily  be  derived  by  combining  the  procedures 
for  (4.62)  and  (4.60). 

A simple  treatment  is  required  for  each  of  the  above  integrals,  in  order 
to  be  able  to  use  a suitable  formula  in  Sec.  3,  after" transformation  from  V 


to  £>d  . Denoting  an  arbitrary  one  of  the  integrals  (4.59)  to  (4.62)  by  T F , 
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we  assume  that  F is  continuous  and  bounded  on  r , and  we  first  set 

(4.63)  F = LF  + E , 

where  LF  is  a simple,  suitable  explicit  interpolation  of  F at  the  end-points 
a and  b of  T . Then  E is  continuous  on  r , and  E(a)  = E(b)  =0  . We 
then  construct  a simple  "polynomial"  p , defined  on  r , such  that 
p(a)  = p(b)  =0  , such  that  P = Tp  can  be  explicitly  expressed,  and  such 
that  the  function  G defined  by 

(4.64)  G = E - p 
satisfies: 

TG  exists  on  T ; 

(4.65)  (TG)  (x)  ->-0  on  x -*■  a or  b , along  r , 

G is  analytic  in  V . 

At  this  point  we  can  apply  a suitable  formula  from  Sec.  3 to  TG  , after 

transformation  from  V to  V,  . 

d 

Let  us  illustrate  the  above  outlined  procedure  on  examples  of  the 
approximation  of  (4.59)  - (4.62). 

EXAMPLE  4.13:  The  Hilbert  Transform  over  a Finite  r . Let  us  make  the 
following  assumptions: 

(i)  F £ B(P)  , where  V is  bounded,  and  T is  finite.  Let  a and  b 
denote  the  end-points  of  T ; 

(ii)  F e Lipa(i?)  , where  0 < a < 1 and  where  V denotes  the  closure 
of  V . 

The  function  LF  referred  to  in  (4.63)  takes  the  form 
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(4.66)  (£F) (x)  = F(a)  + F(b)  . 

The  function  E = F - LF  then  satisfies  <J>'E  <=  B(P)  [36]  , and  moreover,  for 
all  z e V , 

(4.67)  | E (z) | < C| (z-a) (z-b) |a 
where  C is  a constant.  In  particular,  the  integrals 

(4.68)  (#E)  (a)  = |^  Mil  dt  , (flE)(b)  = | dt 

exist  and  are  finite.  Let  p and  q be  polynomials  defined  by 

P(t)  = ^ (t-a)(t-b)(t~|b-ia) 

(b-a)3  3 3 

(4.69) 

Q(t)  - t (t-a)  (t-b)  (t-»a  --w-b  ) 

(b-a)3  3 3 

Then,  define  P and  Q by 

P(x)  = (#p)(x)  = p(x)log(^|)  + (b-a)p'  (x) 

+ j ((b-x)2-(a-x)2}p"(x) 

+ i {(b-x)3-(a-x)3}p"'(x) 

(4.70)  18 

Q(x)  = (»0(x)  - q(x)log(M|)  + (b-a)q'(x) 

• + J ( (b-x)2-(a-x)2}q"(x) 

♦ Tg-  ((b-x)3-(a-x)3}q"'(x) 

The  function  p,  q,  P and  Q have  the  following  properties: 

_____  I 
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i 

P(a)  = q(a)  = p(b)  = q(b)  = 0 
(4.71)  P(a)  = Q(b)  = 1 

P(b)  = Q(a)  = 0 

Now  let  us  define  G by 


(4.72)  G(t)  = E(t)-p(t)(flE)(a)-q(t)flE(b) 

The  function  G has  the  following  properties: 

4>'G€B(P)  ; 

G e Lipa(P)  ; 

(•4  73)  G satisfies  (4.67)  in  V ; 

(ffG)  (x)  exists  for  all  x e ' r 

(tfG)(a)  = ffG(b)  = 0 ; 

N*(G,P)  = sup  N(-^~  , t>)  < * . 

xer  "x 

Upon  replacing  t by  ip  (u)  , we  get 


(4.74) 


>.V.  f (g(i/i(u))  ui*W(u) 

■'R  iKu)-x 


Now  using  (3.40)  on  this  last  integral  yields 


du 

u-<t>(x) 
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(4. 75) 


P.V. 


r 


G(t) 

t-x 


dt  + 2h  l 

k=-oo 


sin2{^-  [<t>(x)-kh]} 

x-z, 

k 


<IM_  [l+e~*d/h] 

2sinh(7rd/h) 


Collecting  the  above  results  and  using  the  (3.49)  approximate  form  of 
(3.47),  we  get  the  following  theorem. 


THEOREM  4.7:  Let  F e B(J)  , where  V is  bounded,  and  let  F € Lip^OtJ)  , 
where  0 < a < 1 . Let  h = [ird/faN)]^"  . Then  for  all  x € r 

IP.V.  I ^-dt  - [F(b)  -F(a)  ] - (LF)  (x)iog(^) 


(4.76) 


-P(x)  (#E)  (a)  - Q(x)07E)(b) 

/ 

N r N G(z.)  ri  , nk-j] 

+ h l { l 1—  [1'C-^-  S(k,h)°4>(x) 

k=-N  \J— N d>'  (z.)  zv-z. 


< Cj_  N^/2e"  (^rdotN) 


1/2 


where  E M F-LF  , LF  is  defined  in  Eq.  (4.66),  G in  (4.72),  p and  q in 
(4.69),  P and  Q in  (4.70)  and  where  depends  only  on  F , d and  a . 


EXAMPLE  4.14:  The  Hilbert  Transform  over  [0,”>1  . Let  us  make  the  following 
assumptions : 

(i)  F e B(P)  , where  V is  defined  in  (4.11); 

(ii)  F e Lip^  {z  : |arg  z|  <d  , 0 < |z|  < p}  , for  some  p > 0 

Let  E and  G be  defined  by 


( 
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E(t)  = F(t)-F(0)/(l+t) 

(4.77) 

G(t)  = E(t) - (tfE) (0)  t /(l+t)z 
The  formula  analogous  to  (4.72)  then  becomes 

2 7T 

f00  crtl  ” kh  kh  sin  (4-  [logx-kh]} 

P.V.  — r— p-  dt  + 2h  l e^Gfe  j 

Jq  t_x  k=-°°  x-e™ 


(4.78) 


< i Jj*.(G?P).  [i+e-ird/h] 

“ 1 sinh(^) 


for  all  x e [0  ,■»]  and  for  all  h > 0 , where 


#00 

G(pe"id  ) 

G(pe  id  ) 

H n L 

• 0 

peld-x 

pe‘ld-x 

ap  > 

(4.79)  N*(G,P)  = sup  ) + i2i£ 

Uo  M peld-x  p( 

Thus  be  choosing  h = [ird/(oN)]d^  , it  follows  that 

p-v-  C ^ dt  - lIs1  i^Sd/x) 


- [ ““t  log(l/x)  + 4-]  (ffE)(0) 

(1+x) L 1 x 


(4.80) 


h l { I e^hG(e-*h)  H*S(k,h)«log 

k— N 'j— N e -ejn  > 


< fjV2e-(TrdaN)' 


for  all  x€  [0,»]  , where  depends  only  on  F,  d and  a . 

We  remark  that  in  applications,  the  integrals  (#E)(a)  and  (ffE) (b)  in 
(4.68)  are  evaluated  by  means  of  formula  (4.32)  and  the  integral 
(tfE)(0)  * [ [E(t)/t]dt  in  (4.80)  is  evaluated  by  means  of  formula  (4.36). 

Jo 


^ .'.****  T,*'^  ‘ +■*** 
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We  remark  also  that  (4.76)  and  (4.80)  provide  a convenient  expression  for 
solving  Hilbert-type  problems  (see  e.g.  [11])  and  integral  equations  with 
Cauchy-type  singularities  over  closed  contours . 


EXAMPLE  4.15:  . The  Integral  U0  ,.F  for  r = [0,1]  . In  this  case  we  want 

EzZ 

to  approximate 

(4.81)  (l/g>yf)(x)  = 


[(x-t)2+y2]'6'!5  f(t)dt 


where  f € B(P)  , V as  in  (4.7),  and  where  F € Lip  (U)  , with  0 < a < 1 . 

r00  (a)k(b).  k a 

Using  the  notation  F(a,b;c;x)  = >,  n x for  the  hvpergeometric 

k"°  (c)k  k! 


wn(x,y) 


function,  it  may  be  shown  that  if  n is  a nonnegative  integer,  then 

1 tn[(x-t)2+y2f0_Js  dt 
0 

2 

- i .n-O.  I 

(4.82) 


xn+1(x2>y2)-6-35F(0+Js,l;S^;-J4  ) 


+ (l-x)n"1[(l-x)2+y2]*e'HF(0^,l;^;-^=4“7  ) 

2 (l-x)z+yz 


Let  us  define  Li,  E and  G by 

(Lf)  (x)  - (l-x)f(O)  + xf  (1) 

(4.83)  E » F - LF 

G(x)  - E(x)  - x(l-x){(ax+b)(l/g^f)(0)  ♦ (cx+d)  (I/g >yf)  (1) } 


where 


. - ;--v  >».  s -» ».  ■ 
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a = (3w2-2w7-2w3-w1) 


-1 


, c = -a 


(4.84)  b * (w3-2w2+w1)(w2-w1)'1a  , d = - (w2-w3)  (w^w^  ‘ia  , 


-1 


wn  = wnf0»y) 

1/2 

By  taking  h = |>d/(aN)]  , we  thus  find  by  proceeding  on  for  (3.51),  that 

j\(x-t)2+y2f3_,5f(t)dt  - w0(x,y)f(0) 


wi (x,y) [f (1) -f (0) ] 


[w2(x,y)-w3(x,y)]  [a(i/g^E)  (0)+c(i/g>yE)  (1)  J 


(4.85) 


[w1(x,y)-w2(x,y)][b(i/g>yE)  (0)+d(tfg>yE) (1)] 


N r N G(z.) 

ill  1 

=-N  ' j = -N  <p'  (z. 


) 


cvyV 

(kh-jh)  2+<J>  (y) 2 


-B-h 


Tj-k(<Ky) , B,h)|‘ 


’S(k,h)o<J)(x) 


< C,  NV2  -(TrdaN) 


1/2 


for  all  x€  [0,1]  , where  <p(x)  - log[x/(l-x)]  , is  defined  in  (2.29)  , 
Zy  in  (4.8)  and  where  depends  only  on  f,  d and  a . In  applications 
(#0  yE)  CO)  and  (£/g  yE)(l)  are  evaluated  by  means  of  formula  (4.32).  It 
may  also  be  convenient  to  approximate  wn(x,y)  by  [Lwn(* ,y) ] (x) 

+ twnfzic»^‘fiwn(,»>r)](z]c^S^k’h-)olog  17x  » where  zk  * ^ tanhCkh/2) 

and  where  [Lwn(* ,y) ] (x)  » (l-x)wn(0,y)+xwn(l,y)  . 
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EXAMPLE  4.16:  The  Integral  (K^f)  for  r = [0,1]  . Let  f satisfy  the 
condition  of  the  previous  example.  We  want  to  approximate 

rl 


(4.86) 


(y)(x)  = 


0 


| x-t | B_1f (t)dt  , 6 > 


0 


As  above, 

(Lf)(x)  = (1-x)  f (0)+xf  (1) 

(4.87)  E = f - Lf 

G(x)  = E(x)  - x(l-x)[(a+bx)(FBE)(0)  + (e+dx)(^E)(l); 


(4.88) 


_ (8*2) (B+3) 
8-1 


. _ -2(8*2) 
8-1 


c * 


(8*1) (8*2) (8*3) 
8-1 


(8*1)  6*2 
6-1 


Then,  by  setting 


yx) 


tn|t-x|3_1dt 


(4.89) 


-2L-  x6+n  ♦ (i-x)s  l (?)  xn'-Ik20!C 
(6+l)n  k-0  k 8+k 


1/2 


for  n * 0,1, 2, 3,  and  taking  h * [ird/foN)]  , we  have 
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) 


I 

I 

l 


8-1 


x-t|M  f(t)dt  - [I0(x)f(0)+{I1(x)-I0Cx)}f(l)] 


(4.90) 


+ [a{I3(x)-I2(x)}+b{I2(x)-I1(x)}] 


+ [c{I3(x)-I2(x))+d{I2(x)-I1(x)}] 


t^'1E(t)dt 


(l-t)8'1E(t)dt 


r(6)cos(-y-)  N ^B-y-0  G(zk) 
" k=-N 


1“6  iif,  ■>  8 

4>  Uk) 


N G(z.) 

+ l 2— 

j*-N  <t>'(z.) 
j^k  j 


zj'zk 

6-1 

t'6cos[(k-j)t]dt}s(k,h)°<j>(x) 

j-k 

J 

!0  > 

< Cx  N3/2e'(irdaN) 


1/2 


for  all  xe  [0,1]  , where  Zj  and  <£(x)  are  defined  in  (4.8),  and  where 
C^is  a constant  depending  only  on  F,  d,  a and  8 . 

EXAMPLE  4.17:  The  Integral  V £ for  T = [0,«]  . Let  f satisfy  the 
condition  in  Ex. 4. 14.  We  shall  approximate  the  integral 

(4.91) 

To  this  end,  we  set 


r 

Jn 


|x-t|6'1f(t)dt 


0 < 8 < 1 


(4.92) 


E a f - Li 


G(x)  - E(x) -a  S-y  (F  E)(0) 

(l+x)Z  8 
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i i 


i- 


where 

(4.93) 

Then,  setting 


„ _ sin(rB) 
a ' ~He — 


x-t 


3 


6-1  dt 
Tn 


BCl+x) 


F(l,B;l+B;— ) + — 1 (1+x)6'1 


1+x  sin(irB) 


(4.94) 


Ix(x)  = f | x-t| 6-1  — — 
JO  (1+t) 


(1+t) - 

_ e+i 

1 x -j  F(2,g;2+B;  - x 2)  + 


6(8+1)  (1+x)"  ’ (1+x)2  sin(TTg)(l+x)2+S 

1/2 

and  taking  h = [7rd/(a'N’>  ] where  o'  * min(a,g)  , we  have 


x-t|S_1f(t)dt  - f(0)I0(x)-a(7gE)(0)I1(x) 


(4.95) 


hr  (6)  cos  (-J)  N , 13  ftkh 

— — l -(4h  - e^Gfe^) 

* k=-N  1 1 B 


N 


+ I e^hG(e^h) 
j*-N 
j’Hc 


Jh  kh 
eJ  -e 

6-1 

J-k 

t Bcos[(k-j)t]dt|’  • 


•S(k,h)olog  x 


< Cj  N^V""*^2 


for  all  x e [0,»]  , where  C.^  depends  only  on  f,  0,  d and  a'  . 

In  applications  (FgE)(0)  is  evaluated  by  means  of  formula  (4.36).  It 
may  also  be  convenient  to  approximate  the  terms  in  IQ(x)  and  Ij(x)  involving 

i 

hypergeometric  functions  by  (4.26).  The  2N  integrals 
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fTT  fl 

t"®cos(kt)dt  = t"^  cos(iTkt)dt  , k=l,2,...  2N  may  be  evaluated 

Jo  J0 

using  (4.32). 


EXAMPLE  4.18:  The  Integral  Q/F)  for  r = [0,1]  ■ 
condition  in  Ex.  4.15.  We  want  to  approximate 


(4.96) 


(VF)(x)  = 


,1 

log|x-t| F(t)dt  . 

0 


Let  F satisfy  the 


We  again  set 


UF)(x)  = (l-x)F(0)+xF(l) 

(4.97)  E = F - LF 

G(x)  = E(x)-(l-x)[(a+bx)(WE)(0)+(c+dx)WE(l)] 


where 


(4.98) 


a = -24 


c = 24 


d = 


78 

5 


Setting 

r1  , n 

I (x)  = log | x-t | t dt 

n J0 

(4.99)  = jij-  {*»+1log  x + (l-x)n+1  log(l-x) 

7n  xk_  \ 

^•k=0  n+l-k  / 

1/2 


for  n 3 0,1,2 ,3,  and  taking  h 3 [ird/(aN)] 


, we  have 


1 
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log|x-t|F(t)dt-[I0(x)f(0)+{I1(x)-I0(x)}f(l)] 


+ [a{I3(x)-I2(x)}+b{I2(x)-I1(x)}]  j E(t) log  t dt 


(4.100) 


+ [c{  1 3 (x)  - 1 2 (x)  }+d{  1 2 (x)  - 1 1 (x) } ] 


E(t)log(l-t)dt 


N f G(z.  ) 

- h I log  1 


k=-N  V(zv)  <J>'(zv) 


* ? i(f£ 

j-N  ♦,Cz.) 
j*=k  -1 


log 


z,  -z. 
,.k  J 


kh-jh 


2n 


k£2s£(k^)tldt 

t 


S(k,h)o$ (x) | 


< Cx  N3/2e'(lTdaN) 


1/2 


where  4>(x)  = log[x/(l-x)]  , and  where  C,  depends  only  on  f,d  and  a 

1 1 rl 


In  applications,  the  integrals 
f1T  1-cos  kt 


E(t)log  t dt  , 


E(t)log(l-t)dt  and 


dt  (k=l,2,. . . ,2N)  are  evaluated  using  (4.32). 


EXAMPLE  4.19:  The  Integral  VF  for  r = f0,°°|  . In  this  case  we  shall  give  an 
approximate  expression  for  the  integral 

* 

.00 

(4. 101)  (J/F)  (x)  * log|x-t|F(t)dt 

Jo 


where  F satisfies  the  conditions  in  Ex.  414.  We  assume  also  that  (VF) (x) 
exists  for  all  x e (0,®)  . Proceeding  as  in  the  previous  examples,  we  set 
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(£F)(x)  = 


(l+t)‘ 


E = F - LF 


(4.102) 


E(x) (a+  -2-  ) 

(1+x)  1+x 


\ [41og  t + 6]E(t)dt 

j0 


h [18-log  t]E(t)dt 
J0 


Let  us  also  set 


Ix(x) 


dt  - ^ log  X 

o (l+tr 


= r tiogix-t 
0 (1+:)3 


(4.103) 


I2(x) 


%X  , . i 1 r-,  . 1 t 

- 777 108  * ^ i*x ll*  a*xK2*x)  > 


0 (1+t) 


1 r,  1 . 0.1  X . 2+X 

5{1'7r^1  log  x T 77^7^  log  i+x 


(1+X)' 


(1+x)- 


. 1 4x+6  ,1^1  (x+2)  -4 


..  -u.  o' 
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1/2 

Then,  by  taking  h = [ud/CaN)]  ' , we  get 

log|x-t|F(t)dt-I0(x)F(0) 

- al^(x)  -bl2(x) 


(4.104) 


kh  e^GCe101) 


N r 

' ' J-H  { 

X ejhG(e'ih)[l08|5ffi4hl 


1 

2¥ 


rir 

0 


l-cost(k-j)t)  dtl|  s(lt>h)olog  x 


<Cj  N3/2e-fe<faN)1/2 


for  all  x s [0,®]  , where  depends  only  on  F,  d and  a 


5.  APPROXIMATION  OF  TRANSFORMS  ON  T . 

The  approximation  of  Fourier  transforms  over  (-“,«)  was  discussed  in 
Sec.  3,  where  the  trapezoidal  rule  yields  the  FFT  method.  In  this  section 
we  briefly  describe  methods  for  approximating  the  semi- infinite  Fourier,  the 
Laplace,  the  Mellin  and  the  Hankel  transforms,  namely 


(5.1) 

,<X> 

F(F,A)  = F(t)elXtdt 

J0 

] 

(5.2) 

L(F,A)  = F(t)e'Xtdt 

Jo 

i 

1 

1 

(5.3) 

Af(F,A)  = f F(t)tX'1dt 

0 

J 

I 

1 

and 

(5.4) 

f°° 

ff  (F,A)  ==  F(t) J (At)dt 

v J0  v 

1 

respectively  [ 24] . 

The  integral  F(F,A)  may  arise  as  a sine  or  cosine  transform,  in  the 

ixt 

process  of  approximating  F(t)e  dt  , when  F has  one  or  more  integrable 

JR 

singularities  on  R . More  generally,  in  the  approximation  of  any  of  the  above 
integrals,  if  F has  a finite  number  of  isolated  singularities  on  the  interval 
of  integration,  then  we  recommend  splitting  up  each  integral  into  a finite 
number  of  integrals,  in  such  a way  that  singularities  occur  only  at  the  end- 
points of  each  interval.  Over  any  such  finite  interval  (a,b)  , we  recommend 
using  the  formula 
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(5.5) 


fb  N _kh 

G(t)dt  2 (b-a)h  j 
U k=-N  (l+e~‘) 


e'~*  rra+be^^ 

wj:  G(t- raH 


1+e 


while  over  (0,®)  we  recommend  using  one  of  the  formulas  (4.36)  or  (4.38). 
Usually  (4.36)  works  well  for  (5.2)  and  (5.3),  while  (4.38)  works  well  for 
(5.1)  and  (5.4). 

An  additional  difficulty  occurs  for  the  application  of  (4.38)  to  the 
approximation  of  (5.1)  and  (5.4),  where  the  integrand  decreases  slowly. 
Consider  for  example,  the  approximation  of 


(5.6) 


by  a trunctation  of 


F(F,A) 


.1/2  sin  At  dt 
t I+t 


1/2 


(5.7) 


S(F,X)  = h l 


k=-«  /1+e'2kh  i+zk 


177  sin(Xzh} 


where  zk  a log[ekh  +/  l+e2kh  ] . This  series  sum  (5.7)  is  quite  accurate 

if  h < min  (1,  m/A)  . For  large  positive  k , -.  kh+log  2 ; hence  if 

A > 1 and  h < m/A  at  least  one  of  the  points  z^  falls  between  every 
consecutive  pair  of  zeros  of  sin  At  , and  S(F,A)  is  then  a very  accurate 
approximation  of  F(F,A)  . We  then  recommend  splitting  the  series  S(F,A) 
into  two  parts  : 

OO  N oo 

h£=h£+h£ 
k»-»  k=-°°  k=N+l 

The  integer  N is  chosen  so  that  \^/ (l+z^  < 1/2  (say)  if  k > N , 
and  so  that  the  series  includes  all  and  only  the  terms  for  which 
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zk  < Po^A  , where  p0  is  a positive  integer.  Direct  summation  and 
truncation  of  the  infinite  series  h ^ ^ Produces  no  difficulty,  since 
this  tail  of  the  series  h ^=oo  converges  very  rapidly.  The  series 
is  summed  by  evaluating  the  first  few 


t = h 
u 


(5.8) 


l 

P*  < (P+1)Tr 

X zk  x 


1/2 


s 


l+e 


TT3T 


1+z, 


sin(Xzk) 


P > Pc 


and  then  applying  Euler's  method  of  summation  to  approximate  the  alternating 

series  7 t 

LP>Po  u 

The  zeros  of  Jy(Xt)  are  asymptotically  equi-spaced  for  large  t and 
the  above  described  procedure  may  be  similarly  applied  to  the  approximation  of 

ffv(F,X)  • 

The  approximation  of  each  of  the  four  integrals  (5.1)- (5.4)  by  the  above 
outlined  procedure  is  described  in  detail  in  [ 24] , where  many  examples  are 
considered,  illustrating  the  accuracy  and  superiority  of  these  methods  over 
other  methods.  It  is  furthermore  shown  in  [24]  that  these  methods  may  be 
used  effectively  for  X up  to  100.  If  X > 100  , we  recomnend  asymptotic 
method  [3,31,41,42,59]. 


6,  APPROXIMATE  SOLUTION  OF  DIFFERENTIAL  EQUATIONS  VIA  THE  SINC-GALERKIN  METHOD. 

In  this  section  we  shall  illustrate  the  application  of  the  cardinal  function 
approximations  to  the  solution  of  second  order  boundary  value  problems  [53]. 

We  defer  the  solution  of  initial  value  problems  to  Sec.  7,  where  we  illustrate 
their  solution  as  a special  case  of  the  solution  of  Volte rra  integral  equations. 

The  method  of  approximate  solution  of  ordinary  and  partial  differential 
equation  boundary  value  problems  which  we  present  here  is  carried  out  by  the 
Galerkin  scheme  [44].  It  is  perhaps  best  illustrated  by  considering  the  solution 
of  the  simple  second  order  linear  boundary  value  problem 

(6.1)  (Lf)(x)  = f"(x)+  u(x)f’(x)+v(x)f(x)-T(x)=0,  xe  r 

(6.2)  f (a)  = f(b)-0  . 

Other  illustrations  of  solutions  of  ordinary  and  partial  differential  equation 
boundary  value  problems  are  given  in  examples  at  the  end  of  this  section. 

Throughout  this  section  we  shall  consider  only  the  solution  of  second  order 
boundary  value  problems. 

Let  A(u)  denote  a diagonal  matrix  with  diagonal  elements  (u _^,u_^+^,. . .iO  , 

T 

i where  u^  = uCz]c)  » 1 denote  the  vector  (1,1,..., 1)  , and  let 

1^  and  1^  denote  the  matrices 
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(6.4)  I(2)  = [6^}] 


Let  us  now  make  some  assumptions  on  Li  in  (6.1).  Using  the  notation 
of  Definition  4.1,  let  y,v  and  a be  analytic  in  V , such  that  (6 . 1) - (6 . 2) 
has  a unique*  solution  f for  which 

(6. 5)  f"/4>'  , f'u/<t>'  , fv/<(>'  and  o/<t>'  e B(P)  . 

and  such  that  f satisfies  (4.15)  on  T . 

We  approximate  f on  V by 

N 

(6.6)  f(x)  a f (X)  = [ f.  S(k,h)°<t>(x)  . 

" k— N K 


ie  assumption  of  uniqueness  may  be  bypassed  via  the  use  of  the 


generalized  inverse. 
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The  Galerkin  scheme  enables  us  to  determine  f^  s f(z^)  by  solving  the  linear 
system  of  equations 


(6.7) 


(£fN,S(k,h)o<t>)  = 0 , k— N,-n+l,. . . ,N  , 


where  the  inner  product  in  (6.7)  is  (nearly  always)  defined  (for  second  order 
problems)  by 


(6.8) 


(u,v)  = 


1 


r <t>’  (x) 


u(x)(  vTx))  dx 


THEOREM  6.1:  Let  the  above  assumptions  be  satisfied.  There  exist  constants 

1/2 

C^,...,Cg  , depending  only  on  f,  d and  a , such  that  if  h - [Ttd/(aN)]  , 
then 


(6.9) 


(6.10) 


(6.11) 


(6.12) 


(6.13) 


(6.14) 


O (x)  ea  t.-v  -a.  /-..■»  j..  u k 


r 4>'(x) 


S(k,h)°<{>(x)dx-h 


< Cj  N-V2e-(«taN) 


1/2 


f *<a.fW  SQc.hWMdx-  H Vk^k> 


r <j>’  (x) 


1/2 


f S(k,h)o<(,Cx)dx-  h.  Uk.f  | < c,  N‘1/2e-(TidaN) 

* » ^ 


r <o’(x) 


uOOf’fo)  S (k ,h)  °<ft (x) dx-  ^ £ f(zj)6^)  | <Cde'(7rdaN) 


1/2 


1/2 


r 4>'(x) 


4>£  j-N 

h 


S(k,h)o<f,(x)dx < C,  N"1/2e‘(7TdaN) 

r*'(x) 

S(k,h)o<(.(x)dx-h  [ f(z.){^  5^  + £ 6^) 

T <*>•  (x)  j-N  3 W*  3K  n 3K  J 


1/2 


< C,  N^-OrdaN) 

" 0 


1/2 
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Let  = (f_N,f_iNJ+1, . . . .f^) 1 be  the  solution  of  the  system 
(6.15)  [J-  I (2)  -A(<t>'-^-)  I (1) +hA(-^)  ] f = hA(-^)l 

A * A * A * 


which  is  obtained  from  (6.7).  Then  the  function  f^  defined  in  (6.6)  then  satisfies 


(6.16) 


| f(x)  ~fN(x)  | < Cj_  N3/2e-(TTdoN)' 


for  all  x e T , where  depends  only  on  f , d and  a . 

We  remark  that  the  approximations  (6.9),  (6.10),  (6.11),  and  (6.13)  may 
be  obtained  by  applying  the  formula  (4.31)  to  the  respective  integrals.  The 
approximations  (6.12)  and  (6.14)  may  be  obtained  by  replacing  f by  f^  and 
then  applying  (4.31).  A different  family  of  approximations  is  also  possible, 
such  as  that  obtained  in  [53].  For  example,  the  above  formulas  yield  the 
approximation 


1 = I ci  S(j  ,h)  °<J>(x)  ]S(k,h)  °<(>(x)dx 

T j J 


* l cj  6jk 


On  the  other  hand,  if  e.g.  f is  bounded  on  T we  find,  after  integration 
by  parts,  and  then  applying  the  above  approximations,  that 

I - - [l  c,S(j,h)o4>(x)]  £ [f(x)S(k,h)»<b(x)]dx 

Jr  J J 

* ] Cj  {j7  f (Xj,4jk)  * 


Both  approximations  have  the  same  order  of  accuracy.  The  approximations  of 
Theorem  6.1  are  usually  simpler  in  form  than  those  we  obtain  after  integration 
by  parts. 


1 
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We  also  remark  that  the  explicit  expressions  (6.9)  - (6.14)  make 

"collocation"  and  "Galerkin"  synonymous  for  this  method. 

The  matrix  Iv  J is  the  dominant  matrix  of  the  system  (6.15).  It  is  a 

symmetric,  negative  definite  matrix,  with  eigenvalues  -X^,k=-N, -N+l, . . . ,N  , 

22  2 (2) 

where  tt  /(N+l)  < < tt  . Thus  J is  a well -conditioned  matrix,  with 

condition  number  less  than  (N+l)^  . The  matrix  1^  is  a skew- symmetric 

matrix  with  determinant  zero.  It  has  eigenvalues  iw^  , where  -tt  < w^  < tt  . 

Contrary  to  the  case  of  finite  difference  or  finite  element  methods  which  lead 

to  sparse  matrices,  the  matrix  in  (6.15)  is  a full  matrix.  However,  the  rate 

of  convergence  (6.16)  of  the  above  method  is  considerably  faster  than  that  of 

finite  difference  or  finite  element  methods,  which  converge  at  the  rate  0(n  ^) 

for  a system  of  order  n , where  for  one -dimensional  problems,  q is  usually  1 

-cn^ 

or  2.  Moreover,  under  the  above  assumptions , the  0(e.  ) rate  of  convergence 

cannot  be  improved  (see  Sec.  9),  regardless  of  the  basis.  Due  to  its  rapid 
convergence,  the  present  method  yields  a desired  accuracy  with  a relatively  small 
system  of  equations. 

The  reduction  in  the  amount  of  work  required  is  considerably  greater  in 
two  and  more  dimensions.  In  application,  the  coefficients  of  differential  equations 
in  p dimensions  are  piecewise  analytic  functions  in  each  variable.  Singular- 
ities of  the  solutions  occur  wherever  the  coefficients  of  the  equations  have 
singularities  but  this  occurs  only  on  p-1  dimensional  surfaces.  Thus  (with 
the  exception  of  inverse  problems,  where  determination  of  the  boundaries  is  more 
difficult)  we  can  determine  a priori  the  points,  or  surfaces  where  the  singular- 
ities occur,  and  using  a system  of  order  n , we  can  achieve  as  approximation 
having  an  Otexpf-cn1^^])  error.  This  should  be  compared  with  methods 
based  on  finite  differences,  or  on  polynomial,  or  finite  element-type  approxi- 
mations, for  which  the  error  is  0(n'C//P)  in  p dimensions. 
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r 

I «■ 

We  do  not  have  a precise  idea  at  this  time,  by  how  much  faster  we  can 
^ solve  partial  differential  equations  by  the  above  method  than  by  finite 

difference  or  finite  element  methods.  The  solutions  of  two-dimensional  "model'' 
problems  in  the  examples  which  follow  seem  to  indicate  that  we  can  get  by  with 
less  than  1/3  of  the  work  required  of  classical  methods  to  achieve  3 places  of 
accuracy,  and  less  than  1/10  of  the  work  it  takes  for  5 places.  Preliminary 
calculations  indicate  the  reduction  in  the  amount  of  work  required  to  solve 
3 and  higher  dimensional  problems  is  considerably  greater,  e.g.  by  a factor  of 
100  in  3 dimensions. 

The  above  method  also  easily  reduces  a nonlinear  equation  to  an  algebraic 
system.  For  example,  for  the  case  of  the  problem 

L(f)  - f'(x)  + G(x,f (x) ,f ' (x) ) = 0 

(6.17) 

f (a)  = f(b)  = 0 , 


if  G(*  ,f  ,f ' )/(b'  e B(P)  , then  we  can  make  the  approximation 


(6.18) 


G(xJfxLill201  S(k,h)o<j)(x)dx 
r 4>’(x) 

. hG(zk,f(zk)  f»(zk))  + o(e.ird/h) 


in  which  we  replace  f (zk)  by  fk  , and  f ' (zk)  by  a linear  combination  of 
the  fk  which  is  given  by  combining  (6.11)  and  (6.12).  The  approximate 
solution  of  (6.7)  then  involves  the  solution  of  a system  of  nonlinear  algebraic 
equations  for  the  fk  , which  is  usually  not  an  easy  problem  to  carry  out. 

Mixed  conditions  at  an  end-point  require  a modification  of  the  form  of 
fN  in  (6.6).  Consider,  for  example,  changing  the  conditions  (6.2)  to 


•'*  ••  -v 
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(6.2)’ 

2 

where  |g|  + 

(see  Sec.  4.) 


f(a)  = a , gf(b)  + y£ ' (b)  = 6 , 

> 0 . In  this  case  (6.6)  is  replaced  by  the  approximation 


(6.6)' 


f(x)  s fN(x)  = c^f  + A(x-a)  ♦ B(x-b)2 


N 

+ l c,(x-b)S(k,h)  °<f>(x) 
k=-N  K 


Clearly,  satisfies  (6.2)'  at  x=a  , since  the  sum  c^(x-b)S(k,h) °<(>(x) 

is  zero  at  x=a  and  at  x=b  . Also,  exists  on  all  of  r , except  at 

x=a  . In  particular,  the  derivative  of  c^(x-b)S(k,h) 04>(x)  is  zero  at 

x=b  . Hence  substituting  f^  for  f in  the  second  equation  in  (6.2)' 
enables  us  to  eliminate  B in  (6.6)',  to  get 


fN(x)  - (x-a)2  + Aw(x) 

N ^‘a  8(b-a)3+2Y(b-a) 

(6.6)" 

N 

+ I c,(x-b)S(k,h)°<j>(x) 

k=-N  K 


where 

(6.6)'"  w(x)  = (x-a)[l-{S(:b'a^}^x'a)  ] 

B(b-a)  +2y (b-a) 


The  expression  (6.6)"  for  f^  involves  2N+2  unknowns:  c.n*c-n+1’ • * * ,CN  and 
A . We  can  thus  carry  out  the  approximate  solution  of  (6.1) -(6. 2)'  by  solving 


the  system  of  2N+2  equations 
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(LfN)  (x)  (x-b)S(k,h)o<t)(x)  — , 

Jr  N 4>’(x) 

s hCLfN)(zk)/0>^2  = o , k=-N , -N+l , . . . ,N 

(6.7)’ 

(If  ) (x)  O&L-  dx 

Jr  4>'  (x) 

-h  l _ p 

k=-N  ^ 


It  is  tempting  to  deal  with  (6.2)’  more  simply  by  means  of  the  approximation 

M 

(6.17)  fN(x)  = u£f  + I ckS(k,h)ob(x) 

k= 


where  u and  the  are  unknown.  Then,  since  the 
and  b ,,  it  is  tempting  to  use  the  approximation 


zk  "bunch  up"  near  a 


(6.18) 


f'(b)  * 


U-a 

b-a 


% 


and  to  take  care  of  the  second  condition  in  (6.2)'  via  this  approximation. 

However,  this  does  not  work.  While  (6.17)  may  be  used  as  an  accurate  approximation 
of  f on  T , the  c^  are  close  to  zero  for  k large  and  positive  (and  also 
for  k large  and  negative).  Hence  although  the  error  | f (z^) -fN(z^) | is 
small. 


(6.19) 


fN(b)-fNUN) 

bzN 


may  be  very  large  (see  Sec.  8.2). 


The  coefficients  in  (6. 9) -(6. 14) , for  various  contour  r corresponding 
to  the  mappings  in  Ex.  4. 1-4. 5 and  the  identity  map  of  Sec.  3 are  given  in 
Table  6.1.  The  entries  are  given  both  as  functions  of  x , into  which  we  may 
substitute  x=Zj,  and  also  as  functions  of  w , into  which  we  may  substitute 
w*kh  . 

TABLE  6.1.  THE  COEFFICIEOTS  IN  EOS.  (6-9) -(6. 13) 


(a) 

[0,1] 

1 X 

l0«r3f 

«-*•  <j;(w)  = *5+^  tanWsw 

(b) 

[-1,1] 

i l+x 

<=■*  = tanh*sw 

(c) 

[0,°°] 

log  X 

•*»  tp(w)  = ew 

(d) 

[0,»] 

log  sinh  x 

— ipM  “ log[ew+^  l+e^w] 

. _ 2w 


2A 

, 4e 

*u'x 


‘2  _ 2w 
x e 


tanh  x 


-sech^x 


x 


0 
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Let  us  illustrate  the  above  method  on  the  approximate  solution  of  some 
"model"  boundary  value  problems.  The  computations  of  these  problems  were  carried 
out  by  Burke  [6]. 


EXAMPLE  6.1: 

(6.20)  e2f"-f+l  = 0 , 0 < x < 1 ; f(0)  = f(l)  = 0 


The  system  (6.15)  for  this  equation  becomes 

(6.22)  [£  l(2)-2I(1)-^-  A(x2(1-x)2)]£  = ~ A(x2(l-x)2)l 

e £ 

Solving  this  system  of  equations  yields  the  approximation 

N 

(6.23)  fN(x)  = l fkS(k,h)o  log  0p^)  ; 

k=-N 

1/2 

Taking  h = .75/N  , N=16  , we  get  a solution  which  is  accurate  to  5 places 

if  £=1/5  and  3 places  if  £=1/10  . Similar  accuracy  could  have  been  obtained 
if  instead  of  solving  (6.20)  we  had  solved 

(6.20)'  e2f’-f+x'1(l-x)'1  = 0 , f (0)  = f(l)  = 0 . 

EXAMPLE  6.2: 

(6.24)  f"  = f-f3/x2  , 0 < x < » ; f (0)  = f(«)  = 0 

This  problem  is  the  radially  symmetric  form  of  the  three  dimensional  nonlinear 
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Klein- Gordon  equation.  Its  solutions  f satisfy  tp’f  e B(V)  , where  V is 
the  region  defined  in  (4.13)  and  where  0(x)  = log  sinh  x . .Moreover,  f is 
bounded  on  [0,®]  , f(x)  = 0(x)  as  x - 0 , f(x)  = 0(e'x)  as  x - <*>  . Hence 
we  may  expect  the  approximation 

N 

(6.25)  f«(x)  = l fv.  S(k,h)o(log  sinh  x) 

k=-N  K 


to  be  accurate.  Substituting  (6.25)  into  the  differential  equations  and  using 
(6.10),  (6.13),  (6.14)  and  the  entries  (d)  of  Table  6.1  we  get  the  nonlinear 
system  of  equations 


2w 


(6.26) 


1+e  1+e 


2w 


1+e^  { log [ew+/(l+e  j ] }" 


where  f = (f_^  f_^+^  , ...,f^  ) . Talcing  h = . 75/N^^  , N=16  and  solving 
this  system  by  Newton's  method  we  are  able  to  get  the  approximation  of  the 
unique  positive  solution  of  the  problem  (6.24)  which  is  accurate  to  5 dec.  on 
l.»]  . The  problem  (6.24)  has  other  solutions,  and  the  system  (6.26) 
i'  rher  solution^  which  approximatee  these. 


• * -1  . (x.y)  e s - (0, 1)*(0, 1) 


j • 0 on  5S 


H • *S/S 


1/2 


Substituting  the  approximation 
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N N 


(6.29)  u^Cx.y)  = l l u.  .S(i,h)o$(x)S(3  ,h)o<j>(y) 

iN  i— N j— N 13 

into  (6.27),  multiplying  by  [<)>' (x)]"1S(k,h)o<()(x)  [(J)' (y)]'1S(£,h)o<j)(y)  , integrating 
over  S and  using  (6.9)  and  (6.13)  yields  the  system 


(6.30) 


where 


(6.31) 


Setting* 

(6.32) 


we  get 
(6.33) 


BU  + UB  = W 


B = - 2hl ('1'* 

U * [ui;j]  , i,j— N,-N+l N 

7 e2w  2w 

w - -h 

(l+e  ) (1+e  ) 

with  E = [e.j]  = [1]  , i , j=-N, -N+l , . . . ,N 


B = T_1AT  , A 


-N 


-N+l 
‘.X 


N 


AY  + YA  = W' 


where 


(6.34) 


-1 


Y = [yy]  = TUT 
W'  = [w!.]  = TWT" 


T 

*From  the  symmetry  of  the  problem,  UB  = BU  and  we  could  therefore  have 
solved  (6.30)  more  directly,  via  the  formula  U = ^B’^V  . However  the  above 
procedure  is  more  general. 
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The  solution  of  (6.33)  is 

(6.35) 

which  yields 

(6.36)  U = T'-’-YtT'1)1 

Substituting  this  result  into  (6.29)  we  get  an  approximate  solution  u^T  which  is 
accurate  to  5 dec.  on  S . Similar  accuracy  obtains  if  the  -1  in  (6.27)  is 
replaced  by  [xfl-x)y(l-y) ] 1 . 

EXAMPLE  6.4: 

(6.37)  uxx  = ut  » € (0,D  *(0,«°) 

(6.38)  u(x,0+)  = sin(irx)  . 

The  exact  solution  of  this  problem  is 

-A 

(6.39)  u(x,t)  = e sin  ttx 
The  Galerkin  approximation 

- 

u^(x,t)  = e simrx 
N 

l u.  .S(i,h)o<f>(x)S(j,h*)°<f>*(t) 
j=-N  1J 

satisfies  the  boundary  conditions  (6.38),  where  * <f>(x)=log[x/(l-x) ] , 4*(t)=log  t . 

*Due  to  the  0(e  ) rate  of  decrease  of  u(x,t)  as  a function  of  t as 
t -*■  « , we  could  have  achieved  greater  accuracy  by  taking  <J>*(t)  = log  sinh  t , 
as  in  Ex.  6.2. 


(6.40) 


N 

l 

i=-N 


w! . 


\.+\. 
i J 


i 
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Substituting  (6.40)  into  (6.37),  we  arrive  at  the  system  of  equations 
(6.41)  BU  + UC  = V 


£h* 

in  which  w is  evaluated  at  kh,  w*  at  Jlh*  , t^  = Z , = %+%  tanh(jh/2)  , 

and  where  h = . 75/N^^  , h*  = .5/NT^  , and  N=16  . The  equation  (6.41)  is 
solved  for  U be  diagonal izat ion  of  B and  C , and  then  proceeding  similarly 
as  in  Ex.  6.3.  The  resulting  approximate  solution  is  accurate  to  4 dec.  on  S . 


7. 


APPROXIMATE  SOLUTION  OF  INTEGRAL  EQUATIONS. 


l 


i 


The  solution  of  linear  integral  equations,  like  the  solution  of  linear 
differential  equations,  is  analytic  in  each  variable  wherever  the  coefficients 
of  the  equation  are  analytic.  Thus  we  can  determine  a priori,  the  regions  on 
which  the  solution  of  a problem  is  analytic.  This  is  usually  the  case  for  most 
nonlinear  integral  equations  arising  in  applications.  It  is  often  more 
difficult  to  determine  the  exact  nature  of  a singularity  and  it  is  in  these 
instances  that  the  methods  of  Sec.  3 and  4 are  particularly  powerful.  In 
this  section  we  illustrate  the  application  of  some  of  the  approximations  of 
Sec.  3 and  4 on  the  solution  of  Volterra  and  Fredholm  integral  equations. 

Basic  to  the  method  of  approximation  is  the  Galerkin  scheme  (see  [18] 
for  a summary  of  this  scheme  for  the  solution  of  Fredholm  integral  equations ; 
the  case  for  more  general  linear  and  nonlinear  equations  is  discussed  in  [44]). 
For  our  purposes,  the  function  S(k,h) °<J)(x)  play  the  most  important  role  in 
this  scheme,  and  for  the  most  important  kernels  arising  in  applications  the 
explicit  approximations  of  Secs.  3 and  4 enable  us  frequently  to  replace  an 
integral  equation  by  a system  of  algebraic  equations  without  performing  any 
numerical  integration.  Such  procedures  have  been  effectively  carried  out  on 
the  numerical  solution  of  one  and  two-dimensional  singular  integral  equations 
in  [ 32,  36]. 

Let  us  consider  the  case  of  a one -dimensional  problem,  such  as 


(7.1) 


f(x)  = (Kf)(x)+g(x)  , X e T 


where  Kf  takes  on  one  of  the  forms 

m (x) 


(7.2) 


£(x,t)f(t)dt 
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or 

(7.3) 


m (x)  = 


£(x,t)f(t)dt 

r 


We  assume  that  P is  a bounded  domain,  that  g is  analytic  in  V , that 
g e Lip^CJJ)  , 0 < a < 1 , and  that  for  any  such  g , (Kg)  has  the  same  properties 
as  g . We  may  then  expect  (7.1)  to  have  a solution  f with  these  properties. 
Assuming  this  to  be  the  case,  we  approximate  f and  g on  T by 

N+l 

f(x)  a f (x)  = l c.  i/>.(x) 

k=-N+l  K K 

(7.4) 

N+l 

gW  3 gN(x)  = l cL  (x) 

N k=-N-l  * K 


where 

^-N-l(x)  = 

(7.5)  il^OO  = S(k,h) °<p(x)  , k=-N,-N+l,. . . ,N 

\ %+l(x)  = (Mf 

1 

We  then  set 

l 

(7.6)  y j (x)  = Qty.)(x)  , j=-N-l,-N,...,N+l 

and  we  approximate  y.  on  T by  (7.4),  namely 

N+l 

C7-7)  UjW  * PjNW  ■ e:k  *kW  ■ 


Upon  substituting  these  approximations  into  (7.1),  we  arrive  at  the  system  of 
equations 
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i 


i 


N+l 

(7.8)  c,  - l e.,c  = d.  , k=-N-l , -N, . . . ,N+1 

K j=-N-l  JK  J K 

for  detemiining  the  . 

The  most  difficult  part  of  the  procedure  is  the  accurate  approximation  of 
Uj (x)  in  Eq.  (7.6).  To  this  end,  the  approximations  in  Secs.  3 and  4 are 

frequently  helpful,  especially  for  the  case  of  important  singular  integral  equa- 

• • 1/2 
tions  arising  in  applications.  Assuming  that  by  taking,  h = [-rrd/(aN)]  ' we 

can  approximate  g in  (7.4)  and  y.  in  (7.6)  to  within  an  error  of 
1/2  1/2 

0(N  exp[- OdaN)  ' ])  , the  resulting  approximation  of  f by  f^  is  accurate 
to  within  an  error  of  OfN^expf-^daN)1^])  . 

Let  us  illustrate  the  solution  of  some  integral  equations  via  examples . 


EXAMPLE  7.1:  A Vol terra  Integral  Equation.  Let  us  consider  the  approximate 
solution  of  the  Volterra  integral  equation. 


(7.9) 


f(x)  = 


[k(t)f (t)+g(t)]dt  + r(x) 

J0 


, X 6 [0,1]  , 


by  use  of  the  formula  (4.58).  Let  V be  defined  by  (4.9),  and  let  k,g  and 
r e B(P)  . We  assume  furthermore  that  r s Lip^fOjl]  , a > 0 , and  that  (7.9) 
has  a solution  f 6 B(0)  n Lipa[0,l]  . The  method  we  shall  describe  for  ob- 
taining an  approximate  solution  of  f in  (7.9)  may  also  be  applied  to  get  an 
approximate  solution  of  the  initial  value  problem 

(7.10)  ^ = A(x)y  + b(x)  , xe  (0,1)  , y(0)  = y0 

3x  " 


in  which  b , y and  yQ  are  vectors  and  A is  a matrix. 
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Let  us  approximate  f and  r on  [0,1]  by 

N+l 

f(x)  a f (X)  = l c.  i|>.(x) 

1 j— N-l  3 J 

N+l 

l Pa  ^(x) 
j— N-l  J J 


r(x)  ^ rN(x) 


(7.11) 


'I'-N-iOO  = 1-x  , %+1(x)  = x 

^j(x)  = S(j,h)°log  ^ , j=-N,-N+l,...,N 

P-N-1  = > PN+1  = 

Pj  = -U-Zj)P.N.1-ZjP  + r(Zj) 


z.  = J5-+J5  tanh  (jh/2)  , 


Substituting  these  expressions  into  (7.9),  using  (4.58)  and  making  the 

r1 

approximation 
equations 


0 


g(x)dx  a h tjm_N  z.  (l-Zj)g(Zj)  , we  get  the  system  of 


C-N-1  = P-N-1 


(7.12) 


N 


N r -I 


,1/2 


Taking  h = [ird/  (aN)  ] and  solving  this  system  for  the  we  get  the 
approximation  f^  in  (7.11),  which  differs  from  f by  less  than  C^N^e 


1/2 
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where  depends  only  on  k,g,r  and  a . 

EXAMPLE  7.2: 

Consider  the  solution  o£  the  problem 

-oo 

(7.13)  f(x)  = k(x-t)f(t)dt+g(x)  , x > 0 , 

J0 

-OO  -CO 

where  ke  HR)  n i}  QR)  and  | g (t)  | “ dt  < <»  | g (t)  | dt  < °°  . Let  us 

Jo  J0 

assume  that  the  transforms 

* . r°° 

(7.14)  K(x)  = elxtk(t)dt  , G.(x)  = elxtg(t)dt 

J]R  j0 

■ 

can  be  explicitly  expressed.  Let  X(x)  * 1 , and  let  d log[l-£(x)]  = 0 . 

J]R 

Then  the  problem  (7.13)  has  a unique  solution  £ on  (0 ,°°)  . For  given  F , 
let  PF  be  defined  as  in  (3.37).  Then  the  Fourier  transform  F+  of  the 
solution  to  (7.13)  may  be  expressed  via  the  formula 

(7.15)  F+  = exp [P4>]  [G++P(G+ [exp{  (1-P)4>}-1] }] 

where  $ = -log[l-K]  . The  formula  (3.48)  may  now  be  used  to  approximate 
P4>  , (1-P)4>  , and  P of  the  remaining  function  in  (7.15).  Using  (3.17),  f(t) 
may  therefore  be  approximated  by  a truncated  Fourier  series  on  (0,  ir/h)  . The 
details  are  carried  out  in  [46,47].  The  convergence  of  this  approximation 
procedure  is  proved  in  [43,46]. 

EXAMPLE  7.3: 

The  approximations  of  this  paper  were  effectively  used  in  [36]  for  ob- 
taining approximate  solutions  of  integral  equations  of  the  form 


-88- 


(7.16) 


f(P)  = X 


^£1  £(Q)dAn  + G(P),P  e S , 
|P-Q|  ^ 


where  S is  a surface  in  ]R  , forming  the  boundary  of  a bounded  region  V . 
For  example,  the  solution  of  the  Neumann  problem  over  a volume  V with  surface 
S 


(7.17) 


Au=0  in  V , |^  = G on  S 


can  be  represented  as  a single  layer  potential 


(7.18) 


u(P)  = 


dA  , p e v . 
S | P-Q|  ™ 


The  unknown  density  function  y then  satisfies  the  integral  equation 

* sr  If  ik  (MMn 


(7.19) 


= S(P) 

2tt 


, P 6 S . 


This  equation  has  eigenvalue  1,  and  from  the  Fredholm  alternative,  we  must  have 

g(Q)dAq  = 0 , in  order  for  a solution  to  exist. 

S 

In  problems  of  the  type  (7.16)  or  (7.19)  arising  in  applications  S is 
a Liapunov  surface,  and  we  generally  expect  to  be  able  to  subdivide  S into  a 
relatively  small  number  of  patches,  S = , so  that: 

(i)  For  each  2,=1,2,...,L  , can  be  parametrized  over  S = [-1,1]  x[-l,l] 
call  the  parametrization  map 


(7.20) 


Tn  ; s * h 
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Let  take  the  form 

(7.21)  T^Cx.y)  = (tp^(x,y)  , t^(x,y)  , t£3^(x,y)) 

where  (x,y)  e S . We  assume  moreover  that  for  each  fixed  y(rsp.x)  e [-1,1]  , 
each  tp-^  (•  ,y)  (rsp.tp-*  (x,°))  e B(P)  n Lip  [-1,1]  , where  a > 0 and  where 
V is  defined  as  in  (4.9); 

(ii)  If  P£  Sj  , then  G(T^)  , satisfies  the  above  conditions  that 
Tp^  satisfies. 

These  assumptions  enable  us  to  reduce  the  equation  (7.16)  over  S to  a 
system  of  integral  equations  over  S . The  approximation  procedure  which  is 
then  applied  is  most  simply  described  by  considering  the  case  L = 1 . In 
this  case  the  resulting  equation  over  S takes  the  form 

(1  ,1  F(x,y;£,n) 

(7.22)  y(x,y)  + y(g,n)  = G(x,y) 

_1  ^/(x-O^fy-n)2 

where  F(»,y;£,ri)  and  G(*,y)  (rsp.  F(x,*;C,n)  and  G(x,0)  belong  to 
B(P)  n Lip  [-1,1]  as  functions  of  x (rsp.  y)  with  the  remaining  variables 
fixed  in  [-1,1]  . 

If  f is  as  arbitrary  continuous  function  defined  on  [-1,1]  , let  us 
set 

N 

(7.23)  L { f(x))  = l f(zk)S(k,h)°q(x) 

k=-N 

where  the  z^  and  <p  are  defined  in  (4.10).  We  then  approximate  y (and 
similarly,  G)  on  S by 
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(7.24) 


where 


N+l 


u(x,y)  a uN(x,y)  = £ y..  <K.(x,y) 

j ,k=-N-l  Jlc 


(7.25) 


4»jk(x,y)  = 

S(j  ,h)  °4>(x) 

S(k 

j,k  = 

-N,-N+l, . . . 

,N 

^-N-l,k(x’y)  = 

(— j^)S(k,h) 

°<Ky) 

Vl,k(x>^  = 

(^)S(k,h) 

°<Ky) 

*j,-N-l(X’y)  = 

(^)S(j  ,h) 

°<f>(x) 

,N+l^x,y^  = 

(l^)S(j,h) 

°<f>(x) 

■N-l,-N-l(‘X,y') 

’-N-l.N+l^’^  = 

(i?)  (i?) 

’N+l.-N-l^’^  = 

%+l,-N-l(x,y^  = 

We  are  thus  led  to  a system  of  (2N+3)2  equation  in  (2N+3)2  unknown. 

1/2 

By  taking  h = [ird/ (aN) ] the  solution  of  this  system  yields  the  approximation 
(7.24)  which  approximates  y on  S to  within  an  error  of  0(Ne '(7TclaN)J5)  . 


'"s  ■"** J " ' 1 • -v‘ 
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For  example,  if  V is  the  unit  ball  and  g(x,y,z)  = a(l-3z  ) on  S , 
where  a = const.,  then  (7.17)  has  the  solution 

(7.26)  u(x,y,z)  = x2+y2-2z2+k 

in  V , where  k is  an  arbitrary  constant.  This  problem  was  solved  in  [ 36  ] 
by  use  of  the  approximation  (7.24),  and  the  resulting  matrix  equation  was 
then  solved  by  use  of  Singular  Value  Decomposition  [15] ; by  taking  N=2  and 
using  the  symmetry  of  the  problem  made  it  possible  to  reduce  the  solution  to 
that  of  a singular  system  of  equations  of  order  13,  yielding  3 places  of 
accuracy. 

EXAMPLE  7.4:  Other  Examples. 

(a)  In  [8]  the  Helmholz  problem  Au  = k u subject  to  Dirichlet  boundary 
conditions  was  solved  on  the  exterior  of  a bounded  region  W in  the  plane, 
via  an  integral  equation  method,  using  the  methods  in  Sec.  4.  It  was  assumed 
that  the  boundary  L of  W consists  of  a finite  number  of  analytic  arcs  , 
with  the  property  that  the  mappings  as  well  as  the  function  g(<jK ) are 

in  B(P)  n Lip  [-1,1]  , where 

(7.27)  <».  : [-1,1]  - L.  , 

where  V is  defined  as  in  (4.9),  and  where  g denotes  tho  boundary  value  of 
u on  L . 

(b)  In  [19]  the  Hilbert  problem 


(7.28) 


F+(t)  - G(t)F_(t)  ♦ H(t)  , t€L 
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was  solved  in  the  complex  plane  via  methods  in  Sec.  4.  It  was  assumed  in  [19] 
that  L consists  of  a finite  number  of  non-overlapping  closed  contours  in  the 
complex  plane,  and  that  these  are  made  up  of  a finite  number  of  analytic  arcs 
Lj  which  can  be  defined  in  the  same  fashion  as  the  Lj  in  (a]  above,  and  such 
that  each  of  the  functions  GC<+>j 3 and  H(<]k)  are  in  B(P)  n Lipa[-l,l]  , where 
V is  defined  in  (4.9). 

(c)  In  [32],  the  problem  of  determining  the  three-dimensional  electric  field 
scattered  by  an  axially  symmetric  body  V in  a plane  wave  was  solved  via  the 
solution  of  an  integral  equation  over  the  surface  S of  B via  the  methods 
of  Sec.  4 of  this  paper.  It  was  assumed  in  [52]  that  the  surface  S is 
described  by 

S = ( (x,y,z)  : x2+y2  = f(z)} 

where  f s B(P)  , V as  in  (4.9),  and  such  that  on  (-1,1) 

2 al  2 a2 

0 < C^l-zl  1 < f(z)  < C2(l-z  ) Z . 

where  and  a2  are  in  (0,1)  . 
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8. 


COMPUTER  IMPLEMENTATION  AND  PITFALLS. 


8.1  Computer  Algorithms. 


Some  of  the  formulas  of  the  previous  section  have  already  been  implemented 
via  computer  programs.  Included  among  these  is  an  automatic  integration  pro- 
gram over  an  arbitrary  interval  (a,b)  using  the  formulas  in  Sec.  4.2  [35], 
programs  for  evaluating  each  of  the  transforms  (5.1)-(5.4)  [22],  programs  for 
the  approximate  solution  of  each  of  the  problems  of  the  examples  in  Sec.  6,  and 
programs  for  the  approximate  solution  of  the  problems  in  Examples  2,  3 and  4 
in  Sec.  7,  as  well  as  in  programs  for  computing  the  solution  to  a Hilbert 
problem. 

Let  us  briefly  describe  the  implementation  of  the  quadrature  program  [35]. 

Let  f 6 B(D)  , and  consider  the  approximation 


(8.1) 


f(x)dx  * T,(f) 

r 


* f (z,  (h)) 

h l — * 

k— « <p'  (zk(h)) 


Let  us  also  set 


(8.2) 


00 


Mu(f)  - h l 
71  k— 


£(z2k-l(h/2)) 

4>’(z2k-i(h/2)) 


9 


so  that 

(8.3)  T(h/2)(£1’I[Th<n®l 

The  error  bound  in  (4.31)  shows  that  when  h is  replaced  by  hh  , the  correct 
number  of  significant  figures  in  the  approximation  (8.1)  double. 
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Assume  then,  that  we  start  with  h=l  (say)  and  then  compute  T^(f)  . Next, 
we  compute  M^(f)  , so  that  if  the  difference 


(8.4) 


Th(f)  - M^f)  = e , 


(8.5)  T(h/2)(f)  = 7(Th(f:i  + Vf)]  = 0(e2)  < e 


In  practice  we  cannot  sum  all  of  the  tenns  in  the  infinite  sums  (8.1)  and 
(8.2).  The  assumption  that 


(8.6) 


_ll2Sl  = o(e‘al(J>W  I)  on  T , a > 0 , 

<t>'  (x) 


then  offers  a convenient  stopping  criteria  in  approximating  the  infinite  sums. 
Suppose  that  we  stop  the  simulations  (8.1)  for  k > 0 when 


(8.7) 


(=  0(6*°^))  < e/3 


Then  we  may  expect  that 


“ f(z^(h))  r -oKh 

h y - < 0(h  7 e oBl) 

k*N+l  <p'  (z^(h))  “ k*N+l 


(8.8) 


. -a(N+l)h 


0(e‘aNh)  - 0(e) 


That  is,  we  may  expect  the  tail  of  the  series  to  be  of  the  same  order  of  magnitude 
as  the  last  included  term.  In  order  to  avoid  stopping  the  algorithm  at  or 
near  a zero  of  f in  practice,  we  make  the  more  reliable  test 


(8.9) 


ifM  , lfW  , lfCW  <e 

4>'(z^)  4> ' C zfj+ 1 ^ ' C ZN+  2) 

A similar  test  is  carried  out  for  negative  k since  the  function  f(x)/<J>'(x) 
may  converge  to  zero  at  different  rates,  as  x+a  or  as  x-*b  along  r . 

These  ideas  form  the  basis  for  the  automatic  integration  algorithm  in  [35]. 
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8.2.  Pitfalls  in  Computation. 


The  accuracy  of  the  formulas  of  this  paper,  in  spite  of  singularities  at  an 
end-point  of  an  interval  is  based  on  our  being  able  to  accurately  compute  the 
function  values  at  the  points  . Consider,  for  example,  the  evaluation  of 


C8.10) 


l> 


x4)  dx 


via  the  formula  (4.34),  in  16  significant  figure  floating  point  arithmetic  [51]. 

The  points  z^(h)  * [e^-l]/[e^+l]  cluster  about  +1  (rsp.-l)  for  k large 

and  positive  (rsp.  negative),  and  the  formula  (4.34)  may  fail,  due  to  roundoff 

2 

error,  resulting  in  the  inaccurate  evaluation  of  [1-z^  (h)]  • For  example, 

if  we  take  h - log2  , k»54  , the  computed  value  of  z^(h)  is 

.99999  99999  99999  0 , so  that  the  computed  value  of  (l-z^2)**5  is 
8 2 -J* 

.707* • • *10  . The  actual  value  of  (1-z^  ) 1 computed  by  means  of  the  formula 
(1-z.2)  **  = (l+e^)/(2e^2)  is  .671*  ••xlO**  . Hence  due  to  roundoff,  the 


(8.11) 


2ekh  u - 2i'h  ■ 2e_kh  fl  _ 2 

VK  T ll'Z.  J * J f 

(l+ej^  k (l+e  "V  k > 


contributes  as  error  of  . 554- • • x10  1 to  the  numerical  approximation. 


(8.12) 


54  -k+l  2k  , 

logh2  l r- j [1"  (~T — ')  1 

k— 54  (l+2*j  2k+1 


of  I . That  is,  it  is  possible  to  achieve  no  more  than  10  significant  figures 
of  accuracy.  If  we  had  carried  the  summation  from  -58  to  58  instead  of 
from  -54  to  54  the  situation  would  have  been  considerably  worse.  In  that 
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58  58 

case  z58  = [2  - 1] / [2  +1]  is  computed  to  be  1.00000  00000  00000  so  that  an 

error  message  results,  since  the  computer  cannot  evaluate  (l.C-1.0)"1*  . 

We  enphasize  that  the  above  difficulty  can  be  easily  remedied,  simply  by 
computing  the  terms  [l-z^2]’*5  by  means  of  the  expression  [l+e^l/^e^2)  . 
Similarly,  the  formula  (see  Eq.  (4.14)) 


(8.13) 


zk  = log[ekh+/  l+e2kh  ] 


is  not  an  accurate  formula  for  computing  z^  for  the  approximations  (4.29) 
and  (4.38)  when  e^1  < .01  , in  that  case  the  formula 


(8.14) 


kh  . 1 3kh  3 Skh  5 7kh 

*k  e 6®  w e • i ne  + ••• 


is  preferable. 

Accurate  computation  of  the  coefficients  is  equally  important  for  the  case 

of  the  formulas  in  Sec.  4.3,  used  for  the  approximation  of  derivatives  over 

finite  and  semi -infinite  intervals.  For  example,  a small  error  in  the  compu- 

2 

tation  of  ffz^)  or  (l-zN  ^ can  cause  a lar8e  error  in  the  approximation 
of  a derivative  in  the  expression 


(8.15) 


N f(zk) 

f(x)  * l *_ 

k— N 2 ' 
(1‘zk) 


(l-xz)  S(k,h)°log 


used  to  approximate  f,f' f'-mj  on  [-1,1]  . 

There  is  one  additional  pitfall  which  we  have  encountered.  For  example, 

let  V be  defined  as  in  (4.9),  let  f € B(P)  , and  let  |f(x)|<  C(l-x2)  on 

% 

(-1,1)  , where  C > 0 . Then  the  approximation 


(8.16) 


f(x)  a l f(zJS(k,h)olog  (ig) 
k—  N K 1 x 
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I 

I 


in  which  h = (md/N)^2  , z ^ = (e^-l)/(e^+l)  is  accurate  for  moderate  values 
of  N if  the  numbers  f (z^)  are  computed  accurately,  and  moreover,  the 
approximation 

-£(zn) 

(8.17)  f'(l)2 


is  then  accurate.  Suppose,  for  example,  that  f is  of  the  order  of  1 , and 
that  the  approximation  (8.16)  is  within  10 of  f for  all  x e [-1,1]  . 

An  error  of  10  ^ in  the  computed  values  of  z^  would  not  change  this  accuracy. 
However,  since  f(x)  -*-0  as  x -*■  ±1  , we  may  have  f(z^)  = ^-*10'6  , and  an 
error  of  10  6 in  f (Zj.)  will  produce  a very  large  error  in  the  approximation 
(8.17).  We  must  therefore  warn  against  using  the  computed  f(z^) — especially 
those  obtained  as  an  approximate  solution  to  a problem  via  the  use  of 
S(k,h)°<J>(x)  as  basis  functions--to  approximate  the  derivatives  of  f at  or 
near  an  end-point  of  an  interval. 


. 
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OPTIMALITY  OF  THE  APPROXIMATIONS. 


The  results  of  this  section  show  that  the  OCe'011  ) rate  of  convergence 
of  the  methods  of  this  paper  cannot  be  improved.  While  the  functions  S(k,h)°4> 
form  a basis  giving  this  rate  of  convergence,  they  are  not  the  only  ones; 
rational  functions  may  also  be  used  to  achieve  this  rate  of  convergence  [48,52] 
The  order  of  the  error  of  approximate  methods  based  on  polynomials  and 
trigonometric  functions  is  well  known  (see  e.g.  [2])  for  many  classes  of 
functions.  We  briefly  cite  some  of  these,  for  purposes  of  comparing  methods  of 
approximation  with  or  without  the  presence  of  a singularity. 

Let  Pn  denote  the  family  of  polynomials  of  degree  < n . 

THEOREM  9.1  [2] : Let  p > 1 , and  let  ep  denote  the  ellipse  with  foci  at 
tl  and  sum  of  semi-ones  equal  to  p . Let  f be  analytic  and  bounded  in  ep 
Then 

(9.1)  inf  sup  | f (x)  -p(x)  | * 0(p*n)  , n-*»  . 

pep  -l<x<l 
n 

That  is,  the  error  of  approximation  converges  to  zero  at  the  0(e'cn) 
rate.  The  rate  of  convergence  (9.1)  is  best  possible  with  regard  to  order. 

The  rate  of  convergence  of  the  error  of  approximation  by  polynomials  is 
considerably  slower  if  a singularity  is  present  at  an  end-point  of  an  interval. 
An  example  of  this  is  illustrated  in  the  following  theorem. 

THEOREM  9.2  [2] : Let  0 < a < 1 . Then  there  exist  positive  constants 


and  C2  such  that 


(9.2) 


“7  < inf  sup  | (l-x2)a-p(x)  | < -J 
n p€Pn  -1<x<1  na 


The  same  drastic  change  in  the  rate  of  convergence  depending  on  whether 
a singularity  is  present  or  absent  occurs  also  for  quadrature  formulas  constructed 
on  the  basis  of  the  formulas  being  exact  for  polynomials  of  a certain  degree. 


For  example,  let  the  numbers  x^ and  , j»l,2,...,n  ; n*l,2,..., 


..  (n)  and  w.fn) 

J J 

be  the  Legendre-Gauss  nodes  and  weights,  so  that  the  approximation 


(9.3) 


f(x)dx  a l w.(n)  f(x.(n)) 
-1  i=l  3 3 


is  exact  wherever  f e p Then  we  have 

Zn-i 


THEOREM  9.3  [40): 


(9.4) 


,1 

-1 


Let  f satisfy  the  conditions  in  Theorem  9.1. 

f (x)dx  - l w.(n)  f(x.  (-n-))  = 0(p'2n)  , n-«°  . 

j-1  3 3 


Then 


The  number  p on  the  right  hand  side  of  (9.4)  cannot  be  replaced  by  a 
smaller  number. 

The  presence  of  a singularity  of  f at  ±1  changes  this  rate  of  convergence 
drastically,  as  illustrated  in  the  following  theorem. 


THEOREM  9.4  flO] : Let  w^ (n)  and  x^ be  defined  as  in  (9.3).  If  a > 0 
and  not  an  integer,  then 


(9.5) 


where  c(o)  depends  only  on  a. 
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For  the  methods  of  this  paper,  the  error  converges  to  zero  at  the  0(e  ) 

rate,  whether  or  not  singularities  are  present  at  an  end-point  of  an  interval. 

. rrfo 

The  best  value  of  c in  this  0(e  ) rate  depends  somewhat,  but  only  mildly, 

on  the  particular  type  of  singularity. 

We  may  therefore  seek  after  methods  which  work  well  in  spite  of  the  presence 
of  a large  class  of  singularities,  and  then  try  to  determine  which,  among  these 
are  the  best. 

Let  us  first  choose  a space  of  functions  with  singularities,  for  purposes  of 
approximation  on  [-1,1]  . 

If  p > 1 , the  space  Hp(U)  is  a convenient  and  well  known  space  of 
functions  in  analytic  function  theory.  Hp(U)  consists  of  the  family  of  all 
functions  f that  are  analytic  in  the  unit  disc  U in  the  complex  plane,  and 


for  which 


(9.6) 


f«  - lim  ^ I2*  |f(rei6)|P  de 

P tm-i  ^ Jn 


Thus  Hp(U)  contains  functions  which  may  or  may  not  have  singularities  at  the 
end-points  of  the  interval  [-1,1]  , such  as  f(x)  - (1-x) ’a(l+x) ’®log(l-x)  , 
where  a ,6  < 1/p  , or  f(x)  - e'x  , etc.  The  closer  p is  to  1 , the 
larger  the  space  Hp(U)  , since,  for  example,  if  l<p'<p  , and  if 
f 6 Hp  (U)  , then  f e , (U)  , whereas  Hp , (U)  has  in  it  functions  that  are 
not  in  Hp(U)  , such  as,  for  example,  the  function  f(x)  - (1-x2)  . 

Next,  given  the  space  Hp(U)  , let  us  consider  the  error  of  approximation 


(9.7) 


1(f)  - <^(f) 


where  f e Hp(U)  , 
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(9.8) 


1(f) 


fl 


f(x)dx  , 0 (f)  - l w (n)f(x.(n)) 
-1  ^ j-l  3 3 


and  where  w.  ^ € I x.^  e U . Let  us  set 
3 3 


(9-9)  % n= 

P’  (n)e(P  (n)CII 
' € (C,x.  e U 


3 


3 


sup 

f€Hp(U),lflp=l 


1(f)  - Q^f) 


The  numbers  ap  n determine  the  best  possible  rate  of  convergence  to  zero  of 

the  quadrature  error.  At  this  time,  the  exact  values  of  the  a and  the 

p,n 

corresponding  quadrature  rules  Qn=Qn*  for  which  | I(f)-Qn*(f)  | < ap  nl|f*p 

for  all  f € Hp(U)  are  not  known.  A number  of  papers  have  been  written  on 

the  estimation  of  upper  bounds  for  a [4,5,17,20,23,50,52,58]  and  some  have 

p ,n 

also  been  written  on  lower  bounds  [4,50,52].  The  results  of  the  following 
theorem  give  the  best  bounds  known  to  date. 


THEOREM  9.5  [52]:  Let  q=p/(p-l)  . Given  any  e > 0 there  exists  an  integer 
n(e)  > 0 such  that  whenever  n > n(e)  , then 

(9.10)  exp[- (51/2Tr+e)n1/2]  < a „ < exp [-{—£-»-  - e}n1/2]  . 

” P’n”  (2qr 

We  remark  that  the  formulas  (4.34)  of  the  present  paper  converge  at  the 
rate  on  the  extreme  right  hand  side  of  (9.10) , for  every  p > 1 , and  the 
formulas  in  [4,48]  also  converge  at  this  rate.  No  formulas  are  known  at  this 
time  which  converge  at  a faster  rate. 

Next,  for  purposes  of  interpolation  on  [-1,1]  , let  p > 1 , let 
Hp*(U)  denote  the  family  of  all  functions  g such  that  f € Hp(U)  , where 
f(z)  - g(z)/(l-z2)  , and  let  Hp*(U)  be  normed  by  lglp  ■ lflp  , where  ■ f I p 
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is  defined  in  (4.6). 
by 


Let  {T  }°°  be  a linear  interpolation  scheme  defined 
n n=l 


(9.11) 


T (g)(x)  = l g(x. 

3-1  J 


(n) 


where  x e u , where  i 

J n,j 

such  that 


is  analytic  in  U for  each  n and 


and 


(9.12)  »Tn(g)«£  ^ C“g“p 

for  all  g e H*fU)  , where  C is  a constant  independent  of  n . Let  us  set 


(9.13) 


inf  sup  sup  | f(x)-T  (f) (x) 

Tn  f e H*(U)  ,1  f«*sl  -l<x<l  n 


THEOREM  9.6  [52] : Given  any  e > 0 , there  exists  an  integer  n(e)>0  , such 
that  whenever  n > n(e)  and  q=p/(p-l)  , 

(9.14)  exp[- (5^^iT+e)n1//2]  < t < exp[-(-V  -e)n1//2] 

p,n  2q^ 

These  upper  and  lower  bounds  are  the  best  ones  known  to  date.  The  formulas 
(4.19)  of  the  present  paper  with  <j>(x)  * log[(l+x)/(l-x)]  converge  at  the 
rate  on  the  extreme  right  of  (9.14).  A rational  function  has  also  been  con- 
structed for  interpolation  of  functions  in  H£  (U)  over  [-1,1]  [52];  this 
also  converges  at  the  rate  on  the  extreme  right  of  (9.14).  No  formulas  are 
known  at  this  time  which  converge  at  a faster  rate. 

Theorem  9. 1-9. 4 show  that  the  formulas  of  this  paper  are  not  as  good  as 
polynomials  in  the  absence  of  singularities,  but  they  are  much  better  when 
singularities  are  present. 
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In  [52]  one  finds  definitions  of  other  Hp  spaces  of  functions  with  the 

property  that  each  of  the  formulas  of  Sec.  4.2  (rsp.  Sec.  4.1)  enjoys  the 

0(exp[-{ — — r -e}  n^])  (rsp.  0(exp[-{-^-  -e}  n^])  rate  of  convergence  in 
(2q)*S  2q^ 

these  spaces. 
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